Principal component pursuit via pcpr
Pattern recognition in environmental mixtures analyses
In this document, we illustrate how to employ Principal Component Pursuit (PCP) for pattern recognition in environmental health (EH) analyses with the pcpr R package. Additional information regarding pcpr’s statistical methodology and computational details are provided in Chillrud et al. (2022).
1 Overview
1.1 The mixtures model
PCP models an observed exposure matrix \(D\) as the sum of three underlying ground-truth matrices: \[ \underset{\text{mixture}}{D_{n \times p}} = \underset{\text{low-rank}}{L_0} + \underset{\text{sparse}}{S_0} + \underset{\text{noise}}{Z_0} \] a low-rank matrix \(L_0\) encoding consistent patterns of exposure, a sparse matrix \(S_0\) isolating unique or outlying exposure events (that cannot be explained by the consistent exposure patterns), and dense noise \(Z_0\). All of these matrices are of dimension \(n \times p\), where \(n\) is the number of observations (e.g., study participants or measurement dates) and \(p\) is the number of exposures (chemical and/or non-chemical stressors). Beyond this mixtures model, the main (EH-relevant) assumption made by PCP is that \(Z_0 \sim \mathcal{N}(\mu, \sigma^2)\) consists of i.i.d. Gaussian noise corrupting each entry of the overall exposure matrix \(D\). A visual example of PCP’s mixtures model is provided below:
Under these minimal assumptions, PCP aims to decompose the observed mixture \(D\) into the above mixtures model via an optimization program, thereby recovering the ground-truth \(L_0\) and \(S_0\) with accurate estimates \(\hat{L}\) and \(\hat{S}\).
1.2 The low-rank matrix
The estimated low-rank matrix \(\hat{L}\) provides information on the consistent exposure patterns, satisfying \[ r = \text{rank}(\hat{L}) \ll \min(n, p). \] The rank \(r\) corresponds to the (relatively few) number of underlying patterns governing the mixture, such as specific sources or behaviors leading to exposure. Notice that \(\hat{L} \in \mathbb{R}^{n \times p}\), meaning it is still defined in terms of the original variables. Put differently, \(\hat{L}\) can be taken as a robust approximation to the true environmental mixture matrix, unperturbed by outliers (captured in \(\hat{S}\)) or noise (handled by PCP’s noise term). In this way, the latent exposure patterns are encoded in \(\hat{L}\) rather than directly estimated. To explicitly obtain the exposure patterns from \(\hat{L}\), PCP may then be paired with various matrix factorization methods (e.g., PCA, factor analysis, or NMF) that yield chemical loadings and individual scores for use in downstream health models.
1.3 The sparse matrix
The estimated sparse matrix \(\hat{S}\) captures unusually high/low outlying exposure events, unexplained by the identified patterns in \(\hat{L}\). Most entries in \(\hat{S}\) are 0, with non-zero entries identifying such extreme exposure activity. The number, location (i.e., support), and value of non-zero entries in \(\hat{S}\) need not be defined; PCP isolates these during optimization.
2 PCP algorithms
root_pcp() |
RRMC() |
|
|---|---|---|
| Convex? | Yes | No |
| Convergence? | Slow | Fast |
| Expected low-rank structure? | Well-defined | Messy |
| Parameters? | \(D, \lambda, \mu\) | \(D, r, \eta\) |
| Supports LOD penalty? | Yes | Yes |
| Supports non-negativity constraint? | Yes | No |
| Rank determination? | Autonomous | User-defined |
| Solver? | ADMM | Iterative rank-based |
2.1 Root PCP
The convex model available in pcpr is Root PCP, denoted here as \(\sqrt{\text{PCP}}\). For a comprehensive technical understanding, we refer readers to Zhang et al. (2021) introducing \(\sqrt{\text{PCP}}\). Here we briefly introduce the root_pcp function in pcpr.
2.2 RRMC
3 Environmental health example analysis: PM2.5 source apportionment
Here we will demonstrate how to apply PCP to environmental mixtures data using an example air pollution source apportionment analysis. Environmental health researchers often aim to identify sources that drive potentially harmful environmental exposures. Investigating potential associations between the identified sources, or patterns of exposure, and adverse health outcomes can then help pave the way toward targeted interventions or public health policy recommendations.
For our analysis, we aim to apportion speciated PM\(_{2.5}\) to its sources using the queens dataset that comes with the pcpr R package. The queens dataset consists of real chemical concentrations (in µg/m\(^3\)) of 26 species of PM\(_{2.5}\) measured every three to six days from 04/04/2001 through 12/30/2021 by an EPA AQS air monitor located in Queens, New York City.
We’ll begin by exploring the raw queens data - corresponding to the \(D\) matrix in our formal mixtures model above - before applying PCP to obtain estimates for \(\hat{L}\) and \(\hat{S}\). We can then compare the PCP-recovered low dimensional structure in \(\hat{L}\) with the observed raw patterns of \(D\).
3.1 Exploring the raw queens data
First, let’s load the pcpr package (along with a few other packages we will need for data processing) and take a look at the queens dataset:
library(dplyr)
library(ggplot2)
library(magrittr)
library(pcpr)
library(stringr)
library(tidyr)
queens ## # A tibble: 2,443 × 27
## Date Al NH4 As Ba Br Cd Ca Cl
## <date> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl>
## 1 2001-04-04 NA 1.62 NA NA NA NA NA NA
## 2 2001-04-07 0 2.66 0 0.012 0.00488 0 0.0401 0.0079
## 3 2001-04-13 0.0094 1.41 0.0016 0.024 0.00211 0.004 0.036 0
## 4 2001-04-19 0.0104 1.22 0.001 0.006 0.00422 0 0.0543 0.003
## 5 2001-04-25 0.0172 0.723 0.0024 0.015 0.00117 0 0.0398 0
## 6 2001-05-01 0.0384 3.48 0.0017 0.041 0.00873 0.001 0.136 0
## 7 2001-05-04 0.0964 6.22 0.0025 0.039 0.0111 0 0.137 0
## 8 2001-05-07 0.004 0.233 0.001 0.016 0.00263 0 0.055 0.0054
## 9 2001-05-10 0.0547 2.04 0.001 0.055 0.00521 0 0.121 0.001
## 10 2001-05-13 0.0215 0.229 0 0.021 0.00122 0 0.0249 0
## # ℹ 2,433 more rows
## # ℹ 18 more variables: Cr <dbl>, Cu <dbl>, EC <dbl>, Fe <dbl>, Pb <dbl>,
## # Mg <dbl>, Mn <dbl>, Ni <dbl>, OC <dbl>, K <dbl>, Se <dbl>, Si <dbl>,
## # Na <dbl>, S <dbl>, Ti <dbl>, NO3 <dbl>, V <dbl>, Zn <dbl>
Let’s spend some time visualizing some of the trends of the raw dataset, to get a sense of what we’re dealing with. We can start by plotting each of the measured chemical species as timeseries. Plotted in black are the raw observed PM\(_{2.5}\) measurements (in µg/m\(^3\)) over time (04/04/2001 - 12/30/2021), and plotted in red are some (rough) lines of best fit:
queens %>%
pivot_longer(
colnames(queens)[-1], names_to = "chem", values_to = "concentration"
) %>%
filter(!is.na(concentration)) %>%
ggplot(aes(x = Date, y = concentration)) +
geom_line() +
geom_smooth(color = "red", formula="y ~ x", method = "loess", span = 0.05) +
facet_wrap(~chem, scales = "free_y") +
labs(x = "Date", y = "Concentration (µg/m^3)") +
theme_bw()Some initial points to make note of, looking at the data:
- The y-axis scale for each chemical’s concentration differs widely (e.g., NH4 ranges from 0 to 9 µg/m\(^3\), while Ti ranges from 0 to 0.04 µg/m\(^3\)). We will have to take this into account during preprocessing (which we do in the next section Preliminary PCA).
- Some chemical species record negative measurements (e.g., Al, Ba, Cd). As explained in the EPA AQS data documentation, this is due to idiosyncrasies in the instruments used to collect the measurements.
- We can see strong seasonal trends in some species (e.g., NO3), as well as a reduction in concentration over time in others (e.g., Ni).
- Many species have prominent outliers, or extreme exposure events (e.g., Cr, Cu). Further, some species have extreme exposure events following seasonal trends, e.g., K.
- Elemental carbon (EC) and organic carbon (OC) are both missing measurements for the first 8 years of the dataset, likely because the air monitors for those two species were not operational from 2001 - 2009. Handling such systematic missingness is out of scope for this tutorial, so we make the simple decision to omit all data from 2001 - 2009, rather than dropping EC and OC, since including both will be helpful in our source apportionment study.
Below we remove the early measurement dates from 2001 to 2009 that are missing observations for EC and OC, yielding the queens_small dataset that we will use moving forward:
start_date_df <- queens %>% select(Date, EC, OC) %>% na.omit() %>% slice_head(n = 1)
start_date <- start_date_df$Date[1]
cat("Start date:", as.character(start_date))## Start date: 2009-04-01
## # A tibble: 1,535 × 27
## Date Al NH4 As Ba Br Cd Ca Cl Cr Cu
## <date> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl>
## 1 2009-04-01 0.044 0.883 0 0 0.0034 0 0.0507 0.677 0 0.0026
## 2 2009-04-04 0.008 0.899 0 0 0.0027 0.009 0.012 0 0.001 0.001
## 3 2009-04-07 0.03 0.623 0 0 0 0 0.0142 0.006 0 0.0039
## 4 2009-04-10 0.016 3.43 0 0 0.0063 0 0.0674 0.0561 0 0.0062
## 5 2009-04-13 0.025 0.774 0.001 0 0.003 0.012 0.0372 0.001 0 0.0025
## 6 2009-04-16 0.002 0.822 0 0.003 0.0051 0.004 0.0458 0 0 0.001
## 7 2009-04-19 0 0.748 0.00105 0 0.0015 0 0.0299 0.0547 0.002 0.001
## 8 2009-04-22 0 2.1 0 0 0.0033 0 0.0488 0.0101 0.002 0.0051
## 9 2009-04-25 0.044 2.47 0 0 0.0037 0 0.0626 0.024 0.001 0.0064
## 10 2009-04-28 0.045 1.33 0 0 0.0043 0 0.0494 0.004 0 0.0023
## # ℹ 1,525 more rows
## # ℹ 16 more variables: EC <dbl>, Fe <dbl>, Pb <dbl>, Mg <dbl>, Mn <dbl>,
## # Ni <dbl>, OC <dbl>, K <dbl>, Se <dbl>, Si <dbl>, Na <dbl>, S <dbl>,
## # Ti <dbl>, NO3 <dbl>, V <dbl>, Zn <dbl>
Next, let’s take a look at the correlation structure of our queens_small data:
queens_small %>%
select(-Date) %>%
as_tibble() %>%
GGally::ggcorr(., method = "pairwise.complete.obs", limits = F, label = F, size = 5)The high dimensionality of our queens_small air pollution mixture matrix gives rise to this relatively complex correlation matrix. No strong patterns jump out right away here, although Na’s correlation with Cl and Mg calls to mind sea salt from the Atlantic Ocean or perhaps road salt kicked up by traffic. We’d like to employ PCP in order to reduce the complexity of our data for more robust downstream analysis. Keep this correlation matrix in mind, since after applying PCP, we’ll examine the correlations matrix of \(\hat{L}\) and compare.
3.2 Preliminary PCA
Before applying PCP to our mixture, it’s good to first take a look at what Principal Component Analysis (PCA) is able to extract from the queens_small data. We do this for a couple of reasons:
- First, we’d like to establish a baseline for comparison. To better understand the effect of applying PCP to the
queens_smalldata, we’d like to see what low dimensional structures exist in the raw data from the outset. - Second, understanding the behavior of the singular values governing the raw data will help inform choices we must make pertaining to PCP’s optimization scheme. A preliminary PCA will enable us to scrutinize our data’s singular values. \(\sqrt{\text{PCP}}\) is well-suited for data exhibiting rapidly decaying singular values (e.g., imaging data), while RRMC is best for more messy data with slowly decaying singular values (environmental mixtures data typically falls into this bucket).
Before we are able to call upon PCA (or PCP, for that matter) we need to first preprocess our data for better numerical stability. In practical terms, both PCA and PCP’s statistical routines are sensitive to data with variable scales. Before arbitrarily normalizing or standardizing our data, let’s examine the distribution of our queens_small dataset:
queens_small %>%
pivot_longer(
colnames(queens_small)[-1], names_to = "chem", values_to = "concentration"
) %>%
filter(!is.na(concentration)) %>%
ggplot(aes(x = concentration)) +
geom_histogram(bins = 50) +
theme_bw() +
facet_wrap(~chem, scales = "free")Most of the chemical species look as though they are normally distributed. We have many choices for how we’d like to preprocess our data, e.g., standardize the data, min-max normalization, etc. Because most of our data appears (roughly) normally distributed, and because in PCA analyses, standardization is typical, let’s go ahead and standardize (scale and center) our data to have \(\mu = 0\), \(\sigma = 1\):
queens_standardized <- queens_small %>%
select(-Date) %>%
scale() %>%
as_tibble()
queens_standardized$Date <- queens_small$Date
queens_standardized %>%
pivot_longer(
colnames(queens_small)[-1], names_to = "chem", values_to = "concentration"
) %>%
filter(!is.na(concentration)) %>%
ggplot(aes(x = concentration)) +
geom_histogram(bins = 50) +
theme_bw() +
facet_wrap(~chem, scales = "free")Great, now we’re ready to run the preliminary PCA analysis! We’ll use the stats::prcomp function for this, which by default omits missing values in our dataset. Instead, let’s handle NAs by imputing with the column-wise means. Below, we define two helper functions, mean_impute and pca:
mean_impute <- function(x) {
x[is.na(x)] <- mean(x, na.rm = T)
x
}
pca <- function(mat, pcs = c("PC1", "PC2", "PC3", "PC4", "PC5", "PC6"), colgroups = NULL, impute = T) {
if ("Date" %in% colnames(mat)) mat <- mat %>% select(-Date)
if (impute) mat <- mat %>% apply(., 2, mean_impute)
#### PCA ####
pca.ln <- prcomp(mat)
#### VARIANCE TABLE ####
singvalues_ln <- matrix(pca.ln$sdev^2)
perc_variance <- round(100*matrix(pca.ln$sdev^2/sum(pca.ln$sdev^2)),1)
pca_summary <- data.frame("Principle component" = 1:min(dim(mat)), "Singular values" = singvalues_ln, "Percent variance" = perc_variance, "Total cumulative variance" = purrr::accumulate(perc_variance, sum))
var_tbl <- kableExtra::kbl(pca_summary, col.names = c("Principle component", "Singular value", "% variance", "Total cumulative variance"), align = "c") %>%
kableExtra::kable_classic(full_width = F, html_font = "Cambria", position = "center") %>%
kableExtra::kable_styling(bootstrap_options = c("hover", "condensed"), fixed_thead = T)
#### LOADINGS ####
pca.ln.ld <- as.data.frame.matrix(pca.ln$rotation)
pca.ln.ld$chem <- row.names(pca.ln.ld)
if (!is.null(colgroups)) {
colgroups <- colgroups %>% dplyr::rename(chem = !!names(colgroups)[1])
} else {
colgroups <- data.frame(chem = colnames(mat), group = "1")
}
grouping <- names(colgroups)[2]
plot_loadings_pca <- pca.ln.ld %>%
tidyr::gather(key = "PC", value = "Loading", -chem) %>%
tibble::as_tibble() %>%
dplyr::right_join(., colgroups, by = "chem")
plot_loadings_pca$chem <- factor(as.character(plot_loadings_pca$chem), levels = unique(as.character(plot_loadings_pca$chem)))
loadings <- plot_loadings_pca %>%
dplyr::filter(PC %in% pcs) %>%
ggplot(aes(x = chem, y = Loading, color = !!sym(grouping))) +
geom_point() +
geom_segment(aes(yend=0, xend = chem)) +
facet_wrap(~ PC) + theme_bw() +
theme(legend.position = "bottom",
axis.text.x = element_text(angle = 45, hjust = 1),
strip.background = element_rect(fill = "white"),
axis.title.x = element_blank(),
axis.title.y = element_blank()) +
geom_hline(yintercept = 0, linewidth = 0.2) +
ggtitle("Principle Component Loadings")
list(var = var_tbl, var_df = pca_summary, load = loadings)
}
prelim_pca <- pca(queens_standardized)Now let’s take a look at our preliminary PCA. We’ll start by examining the singular values and the proportion of variance explained by each component:
| Principle component | Singular value | % variance | Total cumulative variance |
|---|---|---|---|
| 1 | 5.6835249 | 22.3 | 22.3 |
| 2 | 2.3916708 | 9.4 | 31.7 |
| 3 | 2.1466041 | 8.4 | 40.1 |
| 4 | 1.6284138 | 6.4 | 46.5 |
| 5 | 1.4501119 | 5.7 | 52.2 |
| 6 | 1.2056250 | 4.7 | 56.9 |
| 7 | 1.1309203 | 4.4 | 61.3 |
| 8 | 1.0220365 | 4.0 | 65.3 |
| 9 | 0.9368778 | 3.7 | 69.0 |
| 10 | 0.8644932 | 3.4 | 72.4 |
| 11 | 0.7973206 | 3.1 | 75.5 |
| 12 | 0.7896729 | 3.1 | 78.6 |
| 13 | 0.7637169 | 3.0 | 81.6 |
| 14 | 0.6590202 | 2.6 | 84.2 |
| 15 | 0.6328695 | 2.5 | 86.7 |
| 16 | 0.6137555 | 2.4 | 89.1 |
| 17 | 0.4826804 | 1.9 | 91.0 |
| 18 | 0.4510860 | 1.8 | 92.8 |
| 19 | 0.3998378 | 1.6 | 94.4 |
| 20 | 0.3005401 | 1.2 | 95.6 |
| 21 | 0.2833715 | 1.1 | 96.7 |
| 22 | 0.2543039 | 1.0 | 97.7 |
| 23 | 0.2154730 | 0.8 | 98.5 |
| 24 | 0.1757374 | 0.7 | 99.2 |
| 25 | 0.1378358 | 0.5 | 99.7 |
| 26 | 0.0368678 | 0.1 | 99.8 |
prelim_pca$var_df %>%
ggplot(aes(x = Principle.component, y = Singular.values)) +
geom_point() +
geom_segment(aes(y = Singular.values, yend = 0, xend = Principle.component)) +
theme_bw() +
labs(x = "Principal component", y = "Singular value") +
ggtitle("Singular values of preliminary PCA")And next let’s examine the top few PCA loadings, corresponding to each of the first six principal components (PCs), or first six patterns governing our data:
3.3 PCP parameter gridsearch
Now we’re about ready to run our first gridsearch.
We can run a grid search using the default_eta and searching through all ranks 1 through 15 with the following code chunk:
pcp_defaults <- queens_standardized %>% select(-Date) %>% get_pcp_defaults()
default_eta <- pcp_defaults$eta
grid <- grid_search_cv(
mat = D,
pcp_func = RRMC,
grid = data.frame(r = c(15)),
eta = default_eta
)This takes quite a while, so we have saved the results to a .rds file using grid_search_cv’s save_as argument, which we load below:
grid <- readRDS(here::here("my-doc", "huge_eta_search.rds"))
grid$all_stats %>%
group_by(r, eta) %>%
summarize(
avg_rel_err = mean(rel_err),
avg_sparsity = mean(S_sparsity),
avg_rank = mean(L_rank)
) %>%
arrange(avg_rel_err)## # A tibble: 66 × 5
## # Groups: r [6]
## r eta avg_rel_err avg_sparsity avg_rank
## <int> <dbl> <dbl> <dbl> <dbl>
## 1 3 0.2 0.846 100. 3
## 2 3 0.15 0.846 99.9 3
## 3 3 0.1 0.852 99.8 3
## 4 3 0.3 0.853 100. 3
## 5 3 0.25 0.856 100. 3
## 6 2 0.3 0.881 100. 2
## 7 3 0.4 0.882 100. 3
## 8 2 0.25 0.882 100. 2
## 9 3 0.35 0.882 100. 3
## 10 2 0.2 0.884 100. 2
## # ℹ 56 more rows
plot_mat <- function(D) {
heatmaply::heatmaply(D, Rowv=F, Colv=F, showticklabels = FALSE, showtickmarks = FALSE)
}colnames(rrmc_out$L) <- colnames(D)
rrmc_pca <- pca(rrmc_out$L, pcs = paste("PC", 1:3, sep = ""))
rrmc_pca$var| Principle component | Singular value | % variance | Total cumulative variance |
|---|---|---|---|
| 1 | 5.206772 | 55.2 | 55.2 |
| 2 | 2.289758 | 24.3 | 79.5 |
| 3 | 1.942686 | 20.6 | 100.1 |
| 4 | 0.000000 | 0.0 | 100.1 |
| 5 | 0.000000 | 0.0 | 100.1 |
| 6 | 0.000000 | 0.0 | 100.1 |
| 7 | 0.000000 | 0.0 | 100.1 |
| 8 | 0.000000 | 0.0 | 100.1 |
| 9 | 0.000000 | 0.0 | 100.1 |
| 10 | 0.000000 | 0.0 | 100.1 |
| 11 | 0.000000 | 0.0 | 100.1 |
| 12 | 0.000000 | 0.0 | 100.1 |
| 13 | 0.000000 | 0.0 | 100.1 |
| 14 | 0.000000 | 0.0 | 100.1 |
| 15 | 0.000000 | 0.0 | 100.1 |
| 16 | 0.000000 | 0.0 | 100.1 |
| 17 | 0.000000 | 0.0 | 100.1 |
| 18 | 0.000000 | 0.0 | 100.1 |
| 19 | 0.000000 | 0.0 | 100.1 |
| 20 | 0.000000 | 0.0 | 100.1 |
| 21 | 0.000000 | 0.0 | 100.1 |
| 22 | 0.000000 | 0.0 | 100.1 |
| 23 | 0.000000 | 0.0 | 100.1 |
| 24 | 0.000000 | 0.0 | 100.1 |
| 25 | 0.000000 | 0.0 | 100.1 |
| 26 | 0.000000 | 0.0 | 100.1 |
print_patterns <- function(pats, colgroups = NULL, n = 1:6, pat_type = "pat", title = "") {
if (!is.null(colgroups)) {
colgroups <- colgroups %>% dplyr::rename(chem = !!names(colgroups)[1])
} else {
colgroups <- data.frame(chem = rownames(pats), group = "1")
}
grouping <- names(colgroups)[2]
colnames(pats) <- paste0(pat_type, stringr::str_pad(1:ncol(pats), width = 2, pad = "0", side = "left"))
pats.df <- pats %>%
tibble::as_tibble() %>%
dplyr::mutate(chem = colgroups[[1]]) %>%
tidyr::pivot_longer(-chem, names_to = "pattern", values_to = "loading") %>%
dplyr::right_join(., colgroups, by = "chem")
pats.df$chem <- factor(as.character(pats.df$chem), levels = unique(as.character(pats.df$chem)))
loadings <- pats.df %>%
dplyr::filter(pattern %in% paste0(pat_type, stringr::str_pad(n, width = 2, pad = "0", side = "left"))) %>%
ggplot(aes(x = chem, y = loading, color = !!sym(grouping))) +
geom_point() +
geom_segment(aes(yend=0, xend = chem)) +
facet_wrap(~ pattern) + theme_bw() +
theme(legend.position = "bottom",
axis.text.x = element_text(angle = 45, hjust = 1),
strip.background = element_rect(fill = "white"),
axis.title.x = element_blank(),
axis.title.y = element_blank()) +
geom_hline(yintercept = 0, size = 0.2) + ggtitle(title)
loadings
}# for FA
library(psych)
library(GPArotation)
library(ggrepel)
rrmc_fa <- fa(rrmc_out$L, nfactors = 3, n.obs = nrow(rrmc_out$L), rotate = "varimax", scores = "regression")## Warning in cor.smooth(r): Matrix was not positive definite, smoothing was done
## Factor Analysis using method = minres
## Call: fa(r = rrmc_out$L, nfactors = 3, n.obs = nrow(rrmc_out$L), rotate = "varimax",
## scores = "regression")
## Standardized loadings (pattern matrix) based upon correlation matrix
## MR1 MR2 MR3 h2 u2 com
## Al -0.05 0.98 0.20 1 0.00076 1.1
## NH4 1.00 0.09 -0.01 1 0.00036 1.0
## As 1.00 0.04 0.02 1 0.00038 1.0
## Ba -0.18 0.97 -0.14 1 0.00086 1.1
## Br 0.95 0.16 0.27 1 0.00043 1.2
## Cd 0.83 -0.13 -0.55 1 0.00067 1.8
## Ca 0.77 0.60 0.23 1 0.00038 2.1
## Cl 0.10 -0.07 0.99 1 0.00137 1.0
## Cr 0.38 0.92 0.12 1 0.00056 1.4
## Cu 0.69 0.71 -0.16 1 0.00042 2.1
## EC 0.83 0.53 -0.17 1 0.00035 1.8
## Fe 0.56 0.82 -0.06 1 0.00047 1.8
## Pb 0.97 0.04 -0.26 1 0.00042 1.1
## Mg -0.12 0.17 0.98 1 0.00128 1.1
## Mn 0.77 0.64 0.04 1 0.00035 1.9
## Ni 1.00 0.03 0.06 1 0.00039 1.0
## OC 0.77 0.63 -0.14 1 0.00037 2.0
## K 0.90 0.34 0.25 1 0.00038 1.5
## Se 0.97 0.00 -0.22 1 0.00042 1.1
## Si 0.06 0.98 0.16 1 0.00072 1.1
## Na 0.06 0.02 1.00 1 0.00135 1.0
## S 0.95 0.30 0.07 1 0.00032 1.2
## Ti -0.06 1.00 -0.02 1 0.00079 1.0
## NO3 1.00 -0.03 0.05 1 0.00042 1.0
## V 0.99 0.10 0.10 1 0.00038 1.0
## Zn 1.00 0.02 -0.03 1 0.00039 1.0
##
## MR1 MR2 MR3
## SS loadings 14.62 7.63 3.74
## Proportion Var 0.56 0.29 0.14
## Cumulative Var 0.56 0.86 1.00
## Proportion Explained 0.56 0.29 0.14
## Cumulative Proportion 0.56 0.86 1.00
##
## Mean item complexity = 1.3
## Test of the hypothesis that 3 factors are sufficient.
##
## df null model = 325 with the objective function = 523.71 with Chi Square = 798399
## df of the model are 250 and the objective function was 331.99
##
## The root mean square of the residuals (RMSR) is 0
## The df corrected root mean square of the residuals is 0
##
## The harmonic n.obs is 1535 with the empirical chi square 0.1 with prob < 1
## The total n.obs was 1535 with Likelihood Chi Square = 505460.2 with prob < 0
##
## Tucker Lewis Index of factoring reliability = 0.176
## RMSEA index = 1.147 and the 90 % confidence intervals are 1.145 1.15
## BIC = 503626.1
## Fit based upon off diagonal values = 1
## Measures of factor score adequacy
## MR1 MR2 MR3
## Correlation of (regression) scores with factors 1 1 1
## Multiple R square of scores with factors 1 1 1
## Minimum correlation of possible factor scores 1 1 1
loadings <- as.data.frame(cbind(rownames(rrmc_fa$loadings[]), rrmc_fa$loadings[])) %>% rename(Variable = V1)
loadings <- loadings %>% mutate_at(colnames(loadings)[str_starts(colnames(loadings), "MR")], as.numeric)
loadings$Max <- colnames(loadings[, -1])[max.col(loadings[, -1], ties.method = "first")] # should be 2:5
loadings %>% kableExtra::kbl(caption = "Loadings") %>% kableExtra::kable_classic(full_width = F, html_font = "Cambria", position = "center") %>%
kableExtra::kable_styling(bootstrap_options = c("hover", "condensed"), fixed_thead = T) %>% kableExtra::scroll_box(width = "100%", height = "400px")| Variable | MR1 | MR2 | MR3 | Max | |
|---|---|---|---|---|---|
| Al | Al | -0.0516901 | 0.9774711 | 0.2027738 | MR2 |
| NH4 | NH4 | 0.9960585 | 0.0862511 | -0.0080859 | MR1 |
| As | As | 0.9986695 | 0.0438255 | 0.0189229 | MR1 |
| Ba | Ba | -0.1774868 | 0.9735461 | -0.1408622 | MR2 |
| Br | Br | 0.9507836 | 0.1590576 | 0.2651090 | MR1 |
| Cd | Cd | 0.8255029 | -0.1342269 | -0.5475927 | MR1 |
| Ca | Ca | 0.7684500 | 0.5971629 | 0.2291292 | MR1 |
| Cl | Cl | 0.0989384 | -0.0731368 | 0.9917119 | MR3 |
| Cr | Cr | 0.3775964 | 0.9180500 | 0.1185138 | MR2 |
| Cu | Cu | 0.6862436 | 0.7090304 | -0.1610070 | MR2 |
| EC | EC | 0.8327887 | 0.5261587 | -0.1710746 | MR1 |
| Fe | Fe | 0.5617517 | 0.8248882 | -0.0593728 | MR2 |
| Pb | Pb | 0.9652234 | 0.0399142 | -0.2575474 | MR1 |
| Mg | Mg | -0.1198352 | 0.1699795 | 0.9774787 | MR3 |
| Mn | Mn | 0.7711625 | 0.6351324 | 0.0396162 | MR1 |
| Ni | Ni | 0.9978816 | 0.0260774 | 0.0562011 | MR1 |
| OC | OC | 0.7672172 | 0.6266471 | -0.1353504 | MR1 |
| K | K | 0.9043054 | 0.3418836 | 0.2548839 | MR1 |
| Se | Se | 0.9745303 | -0.0042609 | -0.2232732 | MR1 |
| Si | Si | 0.0617839 | 0.9846615 | 0.1609551 | MR2 |
| Na | Na | 0.0639449 | 0.0171711 | 0.9971282 | MR3 |
| S | S | 0.9521291 | 0.2981006 | 0.0653146 | MR1 |
| Ti | Ti | -0.0596583 | 0.9976637 | -0.0178691 | MR2 |
| NO3 | NO3 | 0.9981291 | -0.0343629 | 0.0462689 | MR1 |
| V | V | 0.9901870 | 0.0956123 | 0.1000473 | MR1 |
| Zn | Zn | 0.9992008 | 0.0157969 | -0.0310069 | MR1 |
scores <- as.data.frame(cbind(rownames(rrmc_fa$scores[]), rrmc_fa$scores[])) %>% mutate_all(as.numeric)
scores$Max <- colnames(scores)[max.col(scores, ties.method = "first")]
scores %>%
kableExtra::kbl(caption = "Scores") %>%
kableExtra::kable_classic(full_width = F, html_font = "Cambria", position = "center") %>%
kableExtra::kable_styling(bootstrap_options = c("hover", "condensed"), fixed_thead = T) %>%
kableExtra::scroll_box(width = "100%", height = "400px")| MR1 | MR2 | MR3 | Max |
|---|---|---|---|
| 0.0810764 | -0.2148536 | 4.7956037 | MR3 |
| -0.1758763 | -0.8842593 | -0.4744605 | MR1 |
| -0.5228369 | -0.5805808 | -0.2534992 | MR3 |
| 3.7969090 | -0.1811318 | 0.1175294 | MR1 |
| 0.0170107 | -0.2034000 | -0.3337588 | MR1 |
| 0.3250447 | -0.2863093 | 0.5386281 | MR3 |
| -0.4663021 | 0.0959024 | 0.2860500 | MR3 |
| 2.0423722 | -0.7021320 | -0.2626516 | MR1 |
| 2.2286770 | 0.2290778 | 0.5767719 | MR1 |
| 0.2633727 | 0.6887223 | -0.2459222 | MR2 |
| 1.8155311 | 0.2395212 | 0.7728694 | MR1 |
| 1.0769308 | -0.8702415 | -0.1874851 | MR1 |
| 0.5266173 | 0.3533616 | 0.5098058 | MR1 |
| -0.8160942 | -0.5398569 | -0.4457002 | MR3 |
| 1.2113353 | 0.0682394 | -0.4312160 | MR1 |
| 0.7640474 | -1.2935247 | 0.1144776 | MR1 |
| 0.4466886 | 1.0627672 | -0.1442057 | MR2 |
| 1.0779673 | 1.3788768 | 0.2897110 | MR2 |
| 0.4181035 | 0.9335465 | -0.4760338 | MR2 |
| 0.3750419 | -0.1395150 | -0.5524119 | MR1 |
| 0.2514821 | 0.3204234 | -0.4949604 | MR2 |
| 0.1617738 | 0.0434261 | -0.5974609 | MR1 |
| 0.0420542 | -0.1488882 | -0.6347192 | MR1 |
| 1.2263231 | -0.5701449 | -0.2421441 | MR1 |
| -0.3861805 | -0.6038678 | -0.3965881 | MR1 |
| -0.1971249 | -0.7570192 | 1.8056564 | MR3 |
| -0.7487152 | -0.8230022 | -0.4429393 | MR3 |
| 0.1280251 | 0.1133143 | -0.6027645 | MR1 |
| 0.3094131 | -0.0662733 | -0.5007059 | MR1 |
| 2.0766335 | 0.3415935 | -0.3133643 | MR1 |
| 0.4801154 | 0.3321426 | -0.4596223 | MR1 |
| 1.1938788 | -0.2405081 | 0.0319784 | MR1 |
| -0.3752192 | -0.0473881 | -0.4340776 | MR2 |
| -0.3469356 | -0.5764211 | -0.4363891 | MR1 |
| 0.1316144 | 0.1480555 | -0.3616684 | MR2 |
| 0.7538958 | -0.4999452 | -0.3166927 | MR1 |
| 0.0738393 | -0.6917066 | -0.4745045 | MR1 |
| -0.2114867 | -0.2322134 | -0.5205919 | MR1 |
| 0.9764958 | 0.3764385 | 0.2244427 | MR1 |
| 0.3707445 | 0.2158998 | -0.5639781 | MR1 |
| -0.0026864 | -0.5582212 | -0.1673292 | MR1 |
| 1.7355883 | 0.2183385 | -0.5025919 | MR1 |
| -0.6954272 | -0.4709416 | -0.5881040 | MR2 |
| 0.5312585 | 1.2349601 | -0.5213315 | MR2 |
| 0.9811824 | 0.4687826 | -0.6511636 | MR1 |
| 1.7426769 | 0.1378031 | 0.4326821 | MR1 |
| 1.6559280 | 1.1507623 | -0.0450253 | MR1 |
| -0.0541048 | 1.4125299 | -0.1146895 | MR2 |
| 1.0567550 | 0.4067217 | -0.5909504 | MR1 |
| -1.6517787 | 3.3544503 | -0.1153964 | MR2 |
| -0.7913998 | -0.2904158 | -0.6354265 | MR2 |
| -0.0347403 | 0.1577791 | -0.6542762 | MR2 |
| -0.3987341 | -0.7523689 | 0.0352195 | MR3 |
| -0.4186711 | -0.0036712 | 0.0690681 | MR3 |
| -0.6140767 | -0.5144627 | -0.5385686 | MR2 |
| -0.4940698 | -0.4694065 | 0.0067375 | MR3 |
| -0.8194453 | -0.2367139 | -0.4431056 | MR2 |
| 0.1140868 | -0.1304596 | -0.0188363 | MR1 |
| -0.8874513 | -0.2687159 | -0.5693483 | MR2 |
| 0.4453676 | -0.0974445 | 0.0832252 | MR1 |
| -0.2657420 | -0.4785598 | -0.2827901 | MR1 |
| 1.1043230 | -0.0967694 | -0.3950461 | MR1 |
| -0.4335108 | -0.4511239 | -0.1482413 | MR3 |
| 0.1389915 | -0.3255591 | -0.6686203 | MR1 |
| 0.9104198 | 0.2452368 | -0.5984947 | MR1 |
| -0.6732931 | -0.8213583 | -0.4908518 | MR3 |
| 1.7433415 | 0.7342026 | 0.0122654 | MR1 |
| 3.7598928 | 1.1595188 | -0.4394638 | MR1 |
| -0.5639907 | -0.7303658 | -0.5329985 | MR3 |
| -0.7531444 | -0.8237931 | -0.5765816 | MR3 |
| -0.2864436 | -1.1418485 | 4.6343674 | MR3 |
| 0.5970887 | 0.4595038 | 0.0601781 | MR1 |
| 0.2319998 | -0.6629018 | -0.2543002 | MR1 |
| 6.1124923 | 1.2484408 | -0.1782272 | MR1 |
| -0.8604220 | -0.8825582 | 0.3283628 | MR3 |
| 0.5684062 | 1.2106437 | -0.8168400 | MR2 |
| -0.1544079 | -0.0526669 | -0.1493026 | MR2 |
| 0.9780299 | -0.2324845 | -0.3881576 | MR1 |
| -0.4862931 | -0.9292692 | -0.3154338 | MR3 |
| -0.0614555 | -0.9572409 | -0.5287899 | MR1 |
| 0.3419955 | 1.6183260 | -0.5078970 | MR2 |
| -0.5716393 | -0.9186781 | 0.8424933 | MR3 |
| 1.4938053 | -1.0921639 | -0.1911225 | MR1 |
| 0.7880810 | -0.4325646 | 0.0990783 | MR1 |
| 1.7980323 | -0.0937914 | -0.0899182 | MR1 |
| 2.2078315 | -0.8184787 | -0.3522727 | MR1 |
| 0.7306635 | -0.6848714 | -0.3526078 | MR1 |
| 1.2878815 | -0.6606952 | 0.0515308 | MR1 |
| 0.7837497 | 0.0404739 | -0.3635784 | MR1 |
| 0.2994667 | -0.1671932 | 0.2968346 | MR1 |
| 2.0867267 | -0.6600112 | 0.0539197 | MR1 |
| 0.8232699 | -0.6603968 | 0.1998122 | MR1 |
| 1.2219683 | -0.7524162 | 0.3354554 | MR1 |
| 2.2995283 | -0.6472942 | 0.0994590 | MR1 |
| 2.4249858 | -0.1537253 | 0.3788635 | MR1 |
| 5.6500716 | 0.2592440 | 0.4588566 | MR1 |
| 2.5105456 | 0.2289599 | -0.1965432 | MR1 |
| 2.5272708 | -0.1097335 | -0.4613124 | MR1 |
| 1.0095978 | -0.2681602 | 0.1382686 | MR1 |
| 0.2405660 | -0.6338439 | -0.3383098 | MR1 |
| -0.0925203 | -0.7759681 | -0.0541757 | MR3 |
| 2.1611230 | -0.8010057 | -0.0149140 | MR1 |
| 0.6237412 | -0.3196673 | -0.4936859 | MR1 |
| 0.3157937 | -0.7538374 | -0.0189634 | MR1 |
| 0.4520401 | -0.5387934 | 0.0127517 | MR1 |
| 0.6225843 | -0.5806427 | -0.1626974 | MR1 |
| 1.9551020 | -1.1626441 | -0.3105256 | MR1 |
| -0.0385554 | -0.7986454 | -0.3229689 | MR1 |
| 0.1314752 | -0.1334263 | -0.3948273 | MR1 |
| 0.2460156 | 0.1456481 | -0.1363526 | MR1 |
| 4.2679893 | -1.0302204 | 0.6146390 | MR1 |
| -0.5295391 | -0.6005388 | -0.1809599 | MR3 |
| 0.0088405 | 0.0199748 | -0.3102921 | MR2 |
| 1.2322148 | 0.0297614 | -0.0834964 | MR1 |
| -0.1176026 | -0.8931884 | 1.1135294 | MR3 |
| -0.7400349 | -1.0051236 | -0.5618719 | MR3 |
| 3.1083903 | 2.3694491 | 0.0944544 | MR1 |
| 1.8413000 | 0.3834886 | 0.9146489 | MR1 |
| -0.8362622 | 0.6221428 | -0.3860755 | MR2 |
| -0.3620831 | -0.6678118 | -0.1031946 | MR3 |
| -0.0782141 | -0.1126455 | 0.0760418 | MR3 |
| 1.1430711 | -0.5149046 | -0.1972346 | MR1 |
| 1.5512042 | 1.7262285 | -0.0396571 | MR2 |
| 1.9722680 | 1.8428490 | 3.4201319 | MR3 |
| 0.1824578 | 1.4932779 | -0.1959038 | MR2 |
| 1.1819318 | 1.5910090 | -0.2929092 | MR2 |
| -0.0940834 | -0.9589320 | -0.4221880 | MR1 |
| -0.0045319 | 0.9022000 | -0.4981604 | MR2 |
| -0.6874488 | 2.6864109 | 0.0494737 | MR2 |
| -0.6594410 | -0.3754517 | -0.4839320 | MR2 |
| -0.4304907 | 0.9077372 | -0.2220744 | MR2 |
| 2.0790435 | 1.7950628 | 0.5559775 | MR1 |
| 0.7157535 | 0.8010541 | -0.1879820 | MR2 |
| -0.0389399 | -0.0463574 | 0.0465240 | MR3 |
| -0.0135926 | 0.1161491 | -0.3433647 | MR2 |
| 1.6764088 | -0.3513837 | -0.2661990 | MR1 |
| 0.1706095 | 1.7126882 | -0.9286231 | MR2 |
| 0.1007205 | 0.8664840 | -0.7046681 | MR2 |
| 0.1629021 | -0.6613571 | 0.4145853 | MR3 |
| 0.8336824 | 1.9344255 | 0.2051441 | MR2 |
| -0.1785886 | -0.5792182 | 0.8885868 | MR3 |
| 1.4221290 | -0.0031892 | 0.0522870 | MR1 |
| 0.8570131 | 0.8547118 | -0.5042495 | MR1 |
| -0.4263566 | -0.3332302 | -0.5632275 | MR2 |
| 1.0986987 | -0.3681083 | -0.6706239 | MR1 |
| 0.5873167 | -0.2425405 | -0.1960892 | MR1 |
| 0.1236177 | -0.0664483 | 0.1611528 | MR3 |
| 0.2399607 | -0.2793851 | -0.4348648 | MR1 |
| 0.2964210 | -0.0747946 | -0.3980851 | MR1 |
| -0.2689544 | -0.0567702 | -0.5438300 | MR2 |
| 1.0531220 | 1.0286225 | -0.1111870 | MR1 |
| -0.7045690 | -0.3208099 | -0.4937162 | MR2 |
| 0.9081174 | 1.7188767 | 0.5267676 | MR2 |
| 1.0099941 | 1.2839801 | 0.0013430 | MR2 |
| -0.7786522 | -0.0285106 | -0.1805548 | MR2 |
| 0.1145196 | 0.0346599 | 0.6357661 | MR3 |
| 1.6461233 | 0.1483980 | 0.2888282 | MR1 |
| 0.4845109 | 0.9504312 | -0.3953054 | MR2 |
| -0.6414240 | 0.0151119 | -0.6057159 | MR2 |
| -0.3934751 | 1.2502947 | -0.2557944 | MR2 |
| 0.5104198 | 0.1212605 | 0.1602646 | MR1 |
| -0.5897725 | -0.4526990 | -0.5754374 | MR2 |
| 0.1869337 | 0.1886494 | -0.1074033 | MR2 |
| 1.3712303 | 0.7813277 | -0.6034153 | MR1 |
| -0.5159564 | -0.0352466 | 0.0655931 | MR3 |
| -0.3814615 | -0.9908138 | -0.1205334 | MR3 |
| 0.3627768 | -0.0284462 | -0.6983198 | MR1 |
| -0.4125468 | -0.6916682 | -0.5035996 | MR1 |
| -0.9423026 | -0.9379775 | -0.3117143 | MR3 |
| -0.4885021 | -0.2416682 | -0.7451364 | MR2 |
| 1.3231365 | 1.5250067 | -0.7655793 | MR2 |
| 2.4304159 | 0.8182932 | 0.2399433 | MR1 |
| -0.4418449 | -0.6739767 | -0.4248096 | MR3 |
| -0.2220484 | 0.7661201 | -0.3150759 | MR2 |
| -0.7965088 | -0.2310766 | -0.4055035 | MR2 |
| 1.4152521 | 0.0747046 | -0.7900512 | MR1 |
| -0.0395171 | -0.4326002 | -0.4131414 | MR1 |
| -0.5392898 | -0.6693919 | -0.5630151 | MR1 |
| 0.1878218 | -0.1756450 | -0.1728451 | MR1 |
| -0.6656194 | -0.9971919 | -0.3345037 | MR3 |
| 0.7127787 | 0.2322898 | -0.3580389 | MR1 |
| -0.6879988 | -0.6146689 | -0.6229453 | MR2 |
| -0.6028200 | -0.7276010 | -0.5728823 | MR3 |
| -0.7040649 | 1.6837126 | -0.4578802 | MR2 |
| 1.2543299 | 1.2953016 | -0.4776841 | MR2 |
| 0.6966846 | 1.2413494 | -0.5263223 | MR2 |
| 1.2950087 | -0.3064430 | 4.3561508 | MR3 |
| 2.0583171 | 0.6019274 | -0.5031652 | MR1 |
| 1.4254560 | 0.1922260 | -0.3868243 | MR1 |
| 0.1631737 | -0.5377688 | 0.0011925 | MR1 |
| -0.5628498 | -0.7565792 | -0.2501899 | MR3 |
| -0.1231969 | -0.4363394 | -0.2692232 | MR1 |
| -0.0450966 | -0.5891930 | -0.3805321 | MR1 |
| -0.1794390 | -0.8028956 | -0.3431714 | MR1 |
| 0.1092545 | -0.1054115 | -0.3363054 | MR1 |
| 1.0870164 | 2.6530662 | -0.3362178 | MR2 |
| 1.2032864 | -1.0789124 | 0.9596732 | MR1 |
| 0.6228146 | -0.2529101 | -0.1468190 | MR1 |
| 1.3291072 | -0.2134362 | 0.6583801 | MR1 |
| 0.3845415 | -0.2451617 | -0.3222264 | MR1 |
| 0.2059188 | -0.6486274 | -0.4347729 | MR1 |
| -0.4703339 | -1.1278500 | 3.0577949 | MR3 |
| 0.3396781 | -0.9126818 | -0.2131008 | MR1 |
| -0.2179310 | -0.4411394 | -0.2701635 | MR1 |
| 1.2635961 | -0.0518947 | 0.1141387 | MR1 |
| 0.4979108 | -0.8409490 | -0.2312993 | MR1 |
| 2.5388070 | -0.3565396 | -0.1037327 | MR1 |
| 3.7942105 | -0.7177197 | -0.3821870 | MR1 |
| 0.5860241 | -0.3889473 | -0.5785011 | MR1 |
| 1.0450394 | -0.2936426 | 0.1222715 | MR1 |
| 2.6105817 | -0.8171003 | 0.3236371 | MR1 |
| 8.8558543 | -0.2418781 | 1.8745648 | MR1 |
| 4.5123853 | -0.4452649 | 0.8411076 | MR1 |
| 3.1246607 | -0.6345444 | 0.3664875 | MR1 |
| 0.1319661 | -0.7112776 | -0.2268638 | MR1 |
| 1.3992681 | -0.7919444 | 0.2430627 | MR1 |
| 2.7546988 | -0.4253006 | 1.2059527 | MR1 |
| 3.3019241 | -0.1393201 | -0.4694823 | MR1 |
| 1.5875123 | -0.7800069 | -0.4515860 | MR1 |
| 3.0273585 | -0.3704171 | -0.0180569 | MR1 |
| 0.1682228 | 0.1324543 | -0.0906750 | MR1 |
| 1.5136990 | 0.1558751 | -0.2116527 | MR1 |
| 0.4264055 | -0.0429628 | -0.0200022 | MR1 |
| 2.6171032 | -1.1616082 | 0.1512460 | MR1 |
| 0.7744855 | -0.8086920 | -0.2361216 | MR1 |
| 3.2406949 | -0.0678503 | 0.4788734 | MR1 |
| 0.2678950 | 0.1832925 | 0.0299888 | MR1 |
| 4.3654850 | 1.3684571 | 0.3142868 | MR1 |
| 0.4107358 | -0.5072909 | 1.2979051 | MR3 |
| 1.8644404 | 0.6618966 | 0.7352134 | MR1 |
| 0.7852158 | -0.5567173 | -0.2659279 | MR1 |
| 0.1704379 | 0.0748278 | -0.3201578 | MR1 |
| 0.0744399 | -0.2754846 | 0.4026193 | MR3 |
| -0.2050961 | -0.3373218 | -0.1715237 | MR3 |
| -0.2968660 | -1.3003378 | 8.7136521 | MR3 |
| 0.3844010 | -0.7687465 | -0.0101821 | MR1 |
| 0.2111862 | -0.2692736 | -0.4270621 | MR1 |
| -0.4503988 | -0.6555110 | -0.4449079 | MR3 |
| 0.2598094 | 0.2715178 | -0.3804887 | MR2 |
| 0.3168486 | -0.4956262 | -0.2044324 | MR1 |
| -0.1937094 | -0.4380638 | -0.3182840 | MR1 |
| -0.6849791 | -0.5774941 | -0.4817179 | MR3 |
| 0.2051122 | -0.2665339 | -0.2167265 | MR1 |
| 0.7611227 | -0.1470331 | -0.1448832 | MR1 |
| -0.5894375 | 0.2622500 | -0.5535449 | MR2 |
| -0.4308771 | -0.5824177 | 0.7340618 | MR3 |
| 0.1924050 | 0.7668244 | 0.0459657 | MR2 |
| -0.4465059 | -0.3307522 | -0.5374860 | MR2 |
| 1.8832200 | -0.2010530 | 0.2436684 | MR1 |
| 0.1540543 | -0.8605192 | 0.7404922 | MR3 |
| -0.2975433 | -0.7073806 | -0.4670457 | MR1 |
| 0.8796060 | -0.7108633 | 0.8214977 | MR1 |
| 1.0669886 | 0.0406871 | 0.4411244 | MR1 |
| -0.5371708 | -0.3260840 | -0.4812387 | MR2 |
| -0.5463071 | 0.1868416 | 0.2646801 | MR3 |
| 0.4530652 | -0.4428407 | 1.5779595 | MR3 |
| -0.7583031 | -0.6994451 | 0.6966485 | MR3 |
| 2.5911665 | -0.2758498 | -0.5790147 | MR1 |
| 1.6517177 | 0.2662951 | 0.1043178 | MR1 |
| 2.1542517 | 0.1020534 | 1.4752378 | MR1 |
| 0.4003159 | 0.7852613 | -0.0949435 | MR2 |
| -0.8029629 | 0.4358520 | -0.4223164 | MR2 |
| -0.9014336 | -0.7406472 | -0.1210750 | MR3 |
| 2.5643009 | 1.9701221 | -0.1657794 | MR1 |
| -0.6698272 | -0.3727789 | -0.2103756 | MR3 |
| -0.6412432 | -0.5945272 | -0.2778974 | MR3 |
| 0.2545756 | -0.1921685 | 0.2403679 | MR1 |
| 0.1039545 | 0.0054017 | -0.6418658 | MR1 |
| 1.3321965 | 0.0298186 | -0.4843414 | MR1 |
| 0.1327031 | -0.5323905 | -0.5273875 | MR1 |
| -0.1432602 | 0.4607961 | -0.5181339 | MR2 |
| 0.1462176 | 0.2901067 | 0.1923919 | MR2 |
| 0.1791141 | 0.8381455 | -0.4281444 | MR2 |
| 0.2723080 | 0.2843142 | -0.1663244 | MR2 |
| 0.5973961 | 0.0189579 | -0.0954542 | MR1 |
| 1.3795615 | 0.2042003 | -0.4140796 | MR1 |
| -0.8654603 | -0.2404233 | -0.4456684 | MR2 |
| 0.2587406 | -0.7457814 | 0.2366388 | MR1 |
| 1.8810312 | 0.3784359 | -0.3154466 | MR1 |
| 0.8295934 | 2.0214966 | -0.2378713 | MR2 |
| 0.5933806 | -0.1547266 | -0.2285107 | MR1 |
| 0.9514459 | -0.0965601 | -0.2488947 | MR1 |
| 0.5177101 | 0.7452668 | -0.5030058 | MR2 |
| -0.2349736 | -0.3471821 | -0.1030284 | MR3 |
| -1.0838415 | -0.9058636 | -0.6006952 | MR3 |
| 0.4153132 | 0.0631519 | -0.5668814 | MR1 |
| 0.4513308 | 2.1300343 | -0.3487474 | MR2 |
| -0.4869582 | -0.2293756 | -0.5783945 | MR2 |
| 0.4264398 | -0.0292687 | -0.2443509 | MR1 |
| -0.5064732 | -0.3033510 | -0.5639460 | MR2 |
| -0.1746232 | 0.5307425 | 0.5888776 | MR3 |
| -0.8902565 | -1.0890876 | 2.2615328 | MR3 |
| 0.8510291 | 1.1433297 | -0.2508911 | MR2 |
| 0.2610042 | 0.1287409 | 0.2006784 | MR1 |
| -0.2233585 | 0.3595791 | -0.6315433 | MR2 |
| 1.8797745 | -0.1425529 | 0.2347231 | MR1 |
| 0.9585337 | -0.0606244 | -0.3675955 | MR1 |
| -0.8025637 | -0.5758678 | -0.3142294 | MR3 |
| 1.3160687 | -0.3790952 | 0.0474074 | MR1 |
| 0.3066448 | -0.4399113 | -0.2932203 | MR1 |
| 0.1893144 | -0.2158766 | -0.3885193 | MR1 |
| 0.9881218 | 0.2695058 | 0.3725051 | MR1 |
| -0.3586681 | -0.6773104 | -0.3855770 | MR1 |
| 0.4910487 | 1.0676056 | -0.7568542 | MR2 |
| 2.0310011 | 2.7823966 | -0.5002178 | MR2 |
| -0.7350271 | -0.5610106 | 0.3398527 | MR3 |
| 0.2208306 | -0.7282578 | -0.5468778 | MR1 |
| 1.4725023 | 0.1855094 | -0.3269514 | MR1 |
| -0.2492757 | 0.3071318 | -0.2072599 | MR2 |
| 2.1195293 | -0.5662549 | -0.2666153 | MR1 |
| 0.4280426 | -0.3073045 | -0.4195995 | MR1 |
| 0.3764864 | -0.2814926 | -0.3409435 | MR1 |
| 2.7905255 | 0.2949316 | -0.2441495 | MR1 |
| 0.4806014 | 0.2291202 | -0.3088816 | MR1 |
| 4.3216584 | 2.0017425 | 0.3765289 | MR1 |
| -0.3341741 | -0.7126095 | -0.4711623 | MR1 |
| 0.8272428 | 0.0494874 | 0.6616185 | MR1 |
| -0.2051181 | -0.0793602 | -0.6092156 | MR2 |
| 0.8321643 | -0.4030279 | 0.1137329 | MR1 |
| -0.1500024 | -0.6241267 | -0.2212611 | MR1 |
| 2.6166185 | 0.8912173 | 0.0681703 | MR1 |
| -0.5396299 | -0.8459947 | 2.1160846 | MR3 |
| 1.4885692 | 0.1470437 | -0.3406573 | MR1 |
| 1.9672454 | -0.7197617 | 1.2733181 | MR1 |
| 0.0979584 | -0.8227249 | -0.2669266 | MR1 |
| 0.0160699 | 0.0159664 | 0.0084849 | MR1 |
| 3.0207457 | -0.2176587 | -0.5091983 | MR1 |
| -0.4988283 | -0.3444834 | -0.5413037 | MR2 |
| 0.6935625 | -0.1980137 | -0.1331670 | MR1 |
| 0.3002681 | -0.3272383 | -0.6464131 | MR1 |
| 0.3545758 | -0.5510503 | -0.3844304 | MR1 |
| 0.5009324 | -0.7962695 | -0.2515346 | MR1 |
| 1.3398896 | -0.6383317 | 0.4973805 | MR1 |
| 1.1152154 | -0.2581181 | 0.0338016 | MR1 |
| 2.4105810 | 0.9849416 | 0.1034388 | MR1 |
| 2.4618467 | -0.7247835 | -0.0353963 | MR1 |
| 0.4634358 | -1.0291269 | 0.1703723 | MR1 |
| 0.4919179 | -0.5672552 | -0.0523406 | MR1 |
| -0.3969960 | -0.7901604 | -0.4014775 | MR1 |
| -0.1664245 | -0.9510498 | 0.4151060 | MR3 |
| 0.9504628 | -0.6902031 | -0.4526692 | MR1 |
| 0.0320635 | -0.5577057 | -0.2859690 | MR1 |
| 1.4582289 | -0.0012789 | -0.0756675 | MR1 |
| 0.4374190 | -0.4713628 | -0.4833977 | MR1 |
| 2.1107529 | -0.3354158 | -0.1376034 | MR1 |
| 1.4896503 | -0.7195541 | -0.1388261 | MR1 |
| -0.0588751 | -0.9367446 | 0.1190348 | MR3 |
| 2.8708597 | -0.5981973 | -0.1426819 | MR1 |
| 0.4289547 | -0.8162450 | -0.1887190 | MR1 |
| 0.4694927 | -0.6363940 | 0.2551857 | MR1 |
| 0.3787236 | -0.8884975 | 0.2225618 | MR1 |
| 1.3239086 | -0.1326692 | 0.8506696 | MR1 |
| -0.6130395 | -0.4600221 | -0.4763826 | MR2 |
| 0.0535849 | -0.9388975 | -0.4626166 | MR1 |
| 0.9853654 | -0.8624443 | 2.0306943 | MR3 |
| -0.0096188 | -0.5437191 | -0.3002922 | MR1 |
| 1.4564916 | 0.3295031 | 2.4710578 | MR3 |
| 1.9823658 | 0.1730550 | 0.0667180 | MR1 |
| 2.7117277 | 0.3389181 | 0.7381427 | MR1 |
| 2.2292852 | 0.6008973 | -0.5128264 | MR1 |
| -0.2770764 | -0.3579373 | 1.4260962 | MR3 |
| 0.4266060 | 1.1330990 | -0.0988875 | MR2 |
| -0.6318321 | -0.5741317 | -0.3621861 | MR3 |
| -0.4109593 | 0.5507770 | -0.5711669 | MR2 |
| -0.3917634 | -0.2243658 | -0.3689031 | MR2 |
| -0.7732454 | 0.3624871 | -0.4982371 | MR2 |
| -0.4821846 | -0.4774381 | -0.3576090 | MR3 |
| 1.3296067 | 0.8587096 | -0.1254673 | MR1 |
| -0.3985512 | -0.3925739 | -0.3703114 | MR3 |
| 0.3728600 | -0.6939125 | 0.0668278 | MR1 |
| -0.1761150 | 0.1014128 | -0.2202671 | MR2 |
| -0.5986708 | -0.2664619 | -0.3511008 | MR2 |
| -0.5246346 | -0.2730237 | -0.1453065 | MR3 |
| -0.0246375 | -0.9664023 | -0.1882295 | MR1 |
| -0.3057170 | -0.7174918 | -0.3868334 | MR1 |
| 0.3950703 | 1.3074275 | 0.2433587 | MR2 |
| -0.0045772 | 2.3648379 | -0.3660564 | MR2 |
| -0.2714671 | 0.6701379 | 1.5135811 | MR3 |
| -0.7068120 | -0.1016762 | -0.5721977 | MR2 |
| -0.4811040 | -0.8305230 | 0.0494989 | MR3 |
| 0.1841988 | -1.0152031 | -0.1875763 | MR1 |
| 0.7916770 | 0.7916277 | 0.0004043 | MR1 |
| 0.6440795 | 0.0626273 | -0.5980103 | MR1 |
| -0.3315427 | -0.4702847 | 0.3266129 | MR3 |
| -0.8382346 | -0.6060051 | -0.5651826 | MR3 |
| 0.9737482 | 1.1127593 | -0.5468723 | MR2 |
| -0.2956954 | -0.8863049 | 0.1399787 | MR3 |
| -0.6392716 | -0.4022537 | -0.3479847 | MR3 |
| -0.3975188 | -1.1143751 | 1.2586571 | MR3 |
| 1.2820612 | 0.7159520 | 0.0843342 | MR1 |
| -0.1809196 | 1.2081488 | -0.4978309 | MR2 |
| -0.8724272 | -0.4512083 | -0.5321201 | MR2 |
| 1.4148185 | 1.9465698 | -0.4021645 | MR2 |
| -0.1066475 | 0.1104879 | -0.7184364 | MR2 |
| 0.2053597 | 1.1570003 | -0.0475285 | MR2 |
| -0.2398561 | 0.6934676 | -0.4558340 | MR2 |
| -0.7554921 | -0.3850766 | 0.2401330 | MR3 |
| 0.2217789 | 0.2231727 | -0.4003490 | MR2 |
| 0.6665126 | 0.8833829 | -0.6205105 | MR2 |
| -0.5567030 | -0.5758275 | -0.4762986 | MR3 |
| 0.8355886 | 0.4090435 | 1.2509119 | MR3 |
| 0.5776586 | 0.1737722 | 0.2086669 | MR1 |
| -0.6804386 | -0.8800797 | -0.5938319 | MR3 |
| -0.1298380 | -0.0420613 | -0.2264240 | MR2 |
| 0.0617577 | 0.7418810 | 0.4938471 | MR2 |
| -0.3814735 | -0.4324600 | -0.3730557 | MR3 |
| -0.5862412 | -0.0827423 | 0.7638190 | MR3 |
| 0.7440636 | 0.0348241 | -0.6111988 | MR1 |
| -0.4412016 | -0.1901473 | -0.6506433 | MR2 |
| 0.8898218 | 0.0033916 | -0.4351606 | MR1 |
| 0.2101517 | 0.2164345 | -0.3817779 | MR2 |
| -0.5587222 | 0.9171472 | -0.3218080 | MR2 |
| 1.8050721 | 0.3097683 | -0.3387591 | MR1 |
| -0.7649450 | -0.8959281 | 1.0999218 | MR3 |
| -0.0348666 | 0.4951719 | -0.5674944 | MR2 |
| -0.7516341 | -0.7813485 | -0.5364239 | MR3 |
| 0.7735577 | 0.8454624 | -0.4149363 | MR2 |
| -0.8198114 | -0.1233101 | -0.5892434 | MR2 |
| -0.5028201 | -1.0373416 | 4.4234303 | MR3 |
| -0.7738829 | -0.3363317 | -0.0078759 | MR3 |
| -0.3128359 | 0.1122672 | -0.0668554 | MR2 |
| 1.7148005 | 1.0407961 | -0.1622412 | MR1 |
| -0.6123126 | -0.2861220 | -0.4731721 | MR2 |
| 1.4303771 | -0.4254273 | -0.4745003 | MR1 |
| 0.6500556 | 0.0230060 | -0.4768288 | MR1 |
| -0.8345233 | -0.6296195 | -0.2277709 | MR3 |
| 0.4222216 | 0.6000545 | -0.3012096 | MR2 |
| -0.1629826 | -0.4072309 | 2.2059975 | MR3 |
| 0.2248274 | 0.6620952 | -0.5626615 | MR2 |
| -0.4770124 | -0.1955987 | -0.3972191 | MR2 |
| 0.2078825 | -0.3316803 | -0.5701822 | MR1 |
| -0.4008403 | -0.2791414 | -0.4975582 | MR2 |
| -0.6072510 | -1.5658774 | 4.5790328 | MR3 |
| -0.5115406 | -0.5965736 | -0.6211862 | MR1 |
| -0.3044453 | -0.7110704 | -0.1789408 | MR3 |
| 0.2658161 | -0.1135722 | 0.0894061 | MR1 |
| 2.8654873 | -0.8301279 | -0.1714358 | MR1 |
| -0.1070229 | -0.5939973 | -0.0010695 | MR3 |
| 0.5131546 | -0.1889600 | -0.4417890 | MR1 |
| 0.6852360 | 0.8335433 | -0.5351243 | MR2 |
| 2.5347752 | 0.2066010 | 1.8093561 | MR1 |
| 1.8676804 | 0.1242487 | 0.1638095 | MR1 |
| 3.4003546 | -0.6762088 | 0.0064776 | MR1 |
| 1.9701420 | -0.4571220 | 0.9429971 | MR1 |
| 0.1888883 | -0.6935795 | 0.0533660 | MR1 |
| 3.0986544 | -0.1805636 | 0.3438629 | MR1 |
| 0.3610917 | -0.4575264 | -0.4048670 | MR1 |
| 3.1081957 | 0.5729056 | 0.1160880 | MR1 |
| -0.4643070 | -0.7960470 | -0.5212960 | MR1 |
| -0.1755539 | -0.5705306 | -0.2682120 | MR1 |
| 0.4911862 | -0.6161324 | -0.2281980 | MR1 |
| -0.4265704 | -0.8749363 | -0.1887892 | MR3 |
| 0.4372270 | -0.5296706 | -0.3118105 | MR1 |
| 0.5815561 | -1.0334665 | -0.2559321 | MR1 |
| 1.6208767 | -0.7373584 | -0.3573250 | MR1 |
| 1.6103041 | -0.0973139 | 0.1004733 | MR1 |
| 0.2855423 | -0.1517942 | -0.2667236 | MR1 |
| 0.6135697 | -0.2380171 | -0.2834920 | MR1 |
| 0.8817058 | -0.3259296 | -0.2927349 | MR1 |
| 0.5284210 | -0.7856388 | -0.0352644 | MR1 |
| 0.3926130 | -0.8705455 | -0.1538474 | MR1 |
| 1.8690006 | 0.1090977 | 0.0810668 | MR1 |
| 3.6867721 | -0.9730140 | 1.1198544 | MR1 |
| -0.4033612 | -0.6032385 | 0.9799124 | MR3 |
| 2.0741433 | -1.0877297 | 0.0208275 | MR1 |
| 2.4405135 | -0.3659475 | 0.0372148 | MR1 |
| 0.2537294 | -0.4133149 | -0.0153398 | MR1 |
| 0.3079274 | -0.4598974 | -0.5393489 | MR1 |
| 3.2254295 | -0.8988685 | 0.0272145 | MR1 |
| 0.0610443 | -0.6476733 | 0.5822384 | MR3 |
| 0.4319883 | -0.7922062 | -0.4130892 | MR1 |
| -0.3588245 | -0.6888094 | -0.0846983 | MR3 |
| -0.3326264 | -0.6367731 | -0.5231923 | MR1 |
| 0.3762795 | -0.6113964 | -0.0273926 | MR1 |
| 0.2702245 | -0.5673002 | 0.2024222 | MR1 |
| -0.0616583 | -1.2431581 | 6.9196381 | MR3 |
| -0.4366102 | -0.4619517 | -0.1402601 | MR3 |
| 0.0434018 | -0.6092404 | -0.0863927 | MR1 |
| 0.4488097 | -0.4825546 | -0.3385390 | MR1 |
| -0.6299905 | -0.7225184 | -0.4785131 | MR3 |
| 0.1621407 | -0.0498787 | -0.3210951 | MR1 |
| -0.3954330 | -0.5162193 | -0.5722200 | MR1 |
| 0.6411863 | -0.5972610 | -0.2891695 | MR1 |
| 0.6767521 | -0.3848720 | 0.0735102 | MR1 |
| -0.5850194 | -0.3531472 | 2.5655109 | MR3 |
| 0.6077846 | 0.5838351 | 0.9652006 | MR3 |
| -0.1871280 | -0.6956491 | -0.5291693 | MR1 |
| -0.4090457 | -0.9739595 | 2.9146524 | MR3 |
| -0.6041156 | -0.7024655 | 0.9859122 | MR3 |
| -0.6059308 | -0.0575161 | 1.0397062 | MR3 |
| -0.3820258 | 0.1917835 | 0.0500960 | MR2 |
| 0.3028551 | -0.4040701 | -0.0854046 | MR1 |
| -0.7809564 | 0.1471931 | -0.1978438 | MR2 |
| 0.0072117 | -1.0880731 | 4.5525987 | MR3 |
| -0.6864369 | -0.3808750 | 0.6676885 | MR3 |
| 1.0568603 | 0.5784631 | -0.5305388 | MR1 |
| -0.6926074 | -0.3145410 | -0.4353239 | MR2 |
| 0.5167311 | 2.5559490 | -0.8352747 | MR2 |
| -0.9344095 | -0.8099670 | -0.3522985 | MR3 |
| 1.2292983 | -0.4175363 | 0.9240524 | MR1 |
| -1.3102533 | -0.4491265 | -0.5488724 | MR2 |
| 0.4272781 | -0.0798423 | -0.3407061 | MR1 |
| 0.3024193 | 0.2869051 | -0.1837324 | MR1 |
| 0.9110812 | 0.9590066 | 0.1914096 | MR2 |
| -0.6634591 | 0.1029037 | -0.6321315 | MR2 |
| -0.6649717 | -0.8209099 | -0.3427261 | MR3 |
| -0.0122425 | -0.1314640 | -0.7992740 | MR1 |
| -0.8001158 | -0.4159026 | -0.6692178 | MR2 |
| -0.2163623 | -0.0308981 | -0.2674456 | MR2 |
| 0.2939348 | -0.1567349 | -0.4710862 | MR1 |
| -0.4100732 | 0.2813858 | 0.4342728 | MR3 |
| 0.2953206 | 1.5933362 | -0.5989481 | MR2 |
| 0.3891362 | 0.8919589 | 0.1288275 | MR2 |
| 0.0953011 | 0.0973059 | 2.4179242 | MR3 |
| -1.1489724 | 3.5125120 | 0.8643614 | MR2 |
| -0.5204539 | 0.8882598 | -0.3231638 | MR2 |
| 0.9838370 | 0.3607334 | -0.3278639 | MR1 |
| -0.4706371 | -0.4572573 | -0.5349543 | MR2 |
| -1.5554400 | 7.3598914 | -0.0746122 | MR2 |
| 1.0040178 | 1.1903956 | -0.7159718 | MR2 |
| -0.2249645 | -0.4499434 | -0.5092355 | MR1 |
| -0.2658368 | -0.2245568 | -0.6995156 | MR2 |
| -0.3914531 | 0.3971197 | -0.5855680 | MR2 |
| -0.2236124 | -0.3142222 | -0.6623889 | MR1 |
| 0.0068392 | -0.2478817 | -0.5941856 | MR1 |
| -0.8332191 | -0.4026484 | -0.6639700 | MR2 |
| -0.5519253 | -0.7148858 | 0.6766880 | MR3 |
| -0.8791050 | -0.8138693 | -0.5848689 | MR3 |
| -0.5691432 | -0.8318168 | -0.5762528 | MR1 |
| -0.1263474 | 0.6218829 | -0.5652758 | MR2 |
| 1.1359453 | 0.9253836 | -0.5960043 | MR1 |
| -0.4892888 | 0.0922657 | -0.3156490 | MR2 |
| 0.7096636 | 0.6081041 | -0.5630982 | MR1 |
| 0.2548954 | -0.4262697 | -0.2886507 | MR1 |
| 0.1836689 | -0.9168954 | 1.3003822 | MR3 |
| -0.6483124 | -0.1530236 | -0.7214193 | MR2 |
| -0.3506877 | 0.3886997 | -0.4539800 | MR2 |
| 0.1486302 | -0.0780590 | 0.8686690 | MR3 |
| -0.7125020 | 0.1782676 | -0.6510324 | MR2 |
| -0.4412574 | -0.4248689 | -0.5522029 | MR2 |
| 0.2843531 | 0.8764667 | -0.2164066 | MR2 |
| -0.9038496 | -0.9218189 | -0.5137900 | MR3 |
| -0.4650760 | 1.0562730 | -0.6621910 | MR2 |
| -0.7857417 | -0.1045960 | -0.2747961 | MR2 |
| 0.7037134 | 0.9043301 | -0.7810580 | MR2 |
| 1.4619382 | 0.3694284 | -0.5464058 | MR1 |
| -0.3370479 | -0.8968799 | 0.9601261 | MR3 |
| -0.7622585 | -0.6981956 | -0.0349694 | MR3 |
| -0.6729199 | -0.9498946 | 1.4154259 | MR3 |
| -0.5658091 | -0.4962935 | -0.3627390 | MR3 |
| 0.5034108 | 0.5904840 | -0.8088013 | MR2 |
| -0.2988502 | 0.1249262 | -0.5564089 | MR2 |
| -0.5036095 | -0.2914193 | -0.4687535 | MR2 |
| 0.4454513 | 0.0140003 | -0.4330263 | MR1 |
| 1.0007961 | -0.1444791 | 1.0514077 | MR3 |
| -0.6627267 | -0.8559361 | -0.4686994 | MR3 |
| -0.1191568 | -0.2293769 | 0.2457207 | MR3 |
| -0.6196059 | -0.6215236 | -0.4173956 | MR3 |
| -0.3274522 | -0.4903183 | -0.3977752 | MR1 |
| 2.6155804 | 1.5259915 | -0.3487862 | MR1 |
| -0.6146429 | -0.2105586 | -0.0695593 | MR3 |
| -0.3278302 | -0.5139540 | -0.0477250 | MR3 |
| -0.8603659 | -0.9942669 | -0.4576395 | MR3 |
| -0.3049716 | -0.7492821 | 0.8478632 | MR3 |
| -0.3033708 | -0.4193163 | -0.5090233 | MR1 |
| 4.6006386 | 0.7771618 | -0.5255195 | MR1 |
| 0.1306850 | -0.0642076 | -0.3880180 | MR1 |
| 0.0504284 | -0.6017705 | -0.4695505 | MR1 |
| 0.0126810 | -0.0782172 | -0.2907214 | MR1 |
| 0.2933555 | -1.1585025 | -0.3009115 | MR1 |
| 1.4752463 | -0.8774276 | -0.2341219 | MR1 |
| 0.9318505 | -0.4652353 | 1.2534209 | MR3 |
| -0.1859893 | -0.9029447 | -0.5281640 | MR1 |
| 1.5043671 | 0.2994939 | -0.5714169 | MR1 |
| 0.5097687 | -0.2918122 | -0.1859309 | MR1 |
| 0.1030428 | -0.9968653 | -0.2757833 | MR1 |
| 2.1901565 | -0.6280185 | 0.4729180 | MR1 |
| 2.3107756 | -0.0658734 | 0.5581501 | MR1 |
| -0.0067366 | -0.7659527 | -0.0470963 | MR1 |
| 0.2798988 | -0.1630130 | 0.5590966 | MR3 |
| 1.9557999 | -0.2254804 | 0.2667192 | MR1 |
| 0.5944474 | -0.2783134 | -0.1840601 | MR1 |
| 1.7064477 | -1.0946525 | 0.6535565 | MR1 |
| 0.3094196 | -0.8284977 | 0.0943822 | MR1 |
| 1.6996584 | -0.6852750 | 0.7637515 | MR1 |
| 4.2572231 | 0.1196368 | 0.1496162 | MR1 |
| 1.0836298 | -0.4217711 | -0.0873689 | MR1 |
| 1.9550370 | -0.8264853 | -0.3769512 | MR1 |
| 2.3917282 | -1.0148830 | 0.0555020 | MR1 |
| -0.2329695 | -0.9548698 | 0.1741466 | MR3 |
| 0.2601670 | -0.9011492 | -0.2973881 | MR1 |
| 3.5794779 | -0.5134755 | 1.0781309 | MR1 |
| 0.7867130 | -0.3227058 | -0.2728641 | MR1 |
| 0.7927781 | -0.7106549 | -0.3179059 | MR1 |
| 0.7002427 | 0.1223125 | 0.6574527 | MR1 |
| 0.0663942 | -0.7276819 | 0.0346472 | MR1 |
| -0.0289101 | -0.5224994 | 0.0817423 | MR3 |
| 1.0156439 | -0.8487125 | 0.0226198 | MR1 |
| 2.3465236 | 0.1326478 | 0.9812083 | MR1 |
| -0.1064528 | -0.1776698 | 0.7569763 | MR3 |
| -0.0158771 | 0.0884203 | 0.3520909 | MR3 |
| -0.1061665 | -0.0043838 | -0.2329292 | MR2 |
| -0.3619050 | -0.5566549 | -0.3194921 | MR3 |
| 0.8993796 | -0.0801993 | 0.4015663 | MR1 |
| -0.8087224 | -0.9885950 | -0.4414457 | MR3 |
| 0.9315363 | -0.2980578 | -0.3929936 | MR1 |
| -0.1513723 | -0.9531021 | -0.4324921 | MR1 |
| 0.0020841 | -0.8120692 | -0.2340445 | MR1 |
| 0.5774756 | 1.8628302 | 0.7149784 | MR2 |
| -0.4869962 | -0.6583199 | 7.1900123 | MR3 |
| -0.3162040 | -0.7836640 | 1.4009978 | MR3 |
| -0.4133463 | -0.6323861 | -0.0930567 | MR3 |
| -0.4581217 | -0.0066989 | -0.4456511 | MR2 |
| -0.1719263 | -0.9034481 | -0.4556482 | MR1 |
| -0.7486071 | -1.3078086 | 2.9196634 | MR3 |
| -0.2845979 | -0.0336316 | -0.5866177 | MR2 |
| -0.8772934 | -0.2504583 | -0.5431391 | MR2 |
| 0.9530788 | -0.7908219 | -0.4195736 | MR1 |
| -0.4630977 | 1.2116809 | -0.4554055 | MR2 |
| -0.4296352 | -0.7080944 | 1.6818166 | MR3 |
| -0.3721913 | -0.4312988 | -0.4968577 | MR1 |
| -0.2174026 | 1.6394683 | -0.5529822 | MR2 |
| -0.0321796 | -0.1076343 | -0.5266316 | MR1 |
| 0.4774448 | 0.7109474 | -0.4918060 | MR2 |
| -0.7293523 | -0.9083715 | 0.8602859 | MR3 |
| -0.4474068 | -0.4783802 | 0.1603149 | MR3 |
| 0.8225631 | -0.2227522 | -0.1232435 | MR1 |
| -0.8222980 | 0.1582538 | -0.7037769 | MR2 |
| 0.2583197 | -0.8535695 | 0.5499135 | MR3 |
| -0.4046263 | -0.5892365 | -0.2578605 | MR3 |
| -0.1765754 | -0.5273747 | -0.5281735 | MR1 |
| -0.3808047 | 0.3541756 | -0.5847578 | MR2 |
| -0.2636800 | -0.1383793 | -0.4521848 | MR2 |
| -0.4072722 | -0.3673439 | 0.6480002 | MR3 |
| -0.5367209 | -0.2306676 | -0.6130465 | MR2 |
| -0.4894894 | -0.3817725 | 1.1351655 | MR3 |
| -0.6609041 | 0.4669112 | 0.4717062 | MR3 |
| -0.1686917 | 1.3100344 | -0.3297649 | MR2 |
| -0.4369499 | 0.1989784 | -0.6652121 | MR2 |
| -0.5580363 | -0.2527273 | 0.8399260 | MR3 |
| -0.0660635 | 0.1139718 | -0.6323886 | MR2 |
| -0.2863153 | 0.0814454 | -0.4537512 | MR2 |
| 0.1083223 | 0.0743975 | -0.2376906 | MR1 |
| -0.6489518 | -0.1718607 | -0.4846167 | MR2 |
| -0.7096420 | -0.6498248 | -0.3780756 | MR3 |
| 0.3233675 | 1.1787026 | 0.0819058 | MR2 |
| 0.0506680 | -0.8782312 | -0.4303504 | MR1 |
| -0.0438525 | -0.1858231 | -0.6681106 | MR1 |
| -0.4339301 | -0.0707561 | -0.6120416 | MR2 |
| -0.7042274 | -0.9164574 | 0.0730948 | MR3 |
| -0.9341437 | -0.1306519 | -0.6524116 | MR2 |
| -0.5295593 | -0.0540278 | -0.6398284 | MR2 |
| -0.4306774 | -0.6876896 | -0.2923046 | MR3 |
| -0.8127369 | -0.3263283 | -0.1568014 | MR3 |
| 0.2175071 | 0.1684323 | -0.4287906 | MR1 |
| -0.3454190 | -0.2403823 | -0.4154220 | MR2 |
| -0.8252862 | 4.8530137 | -0.1028609 | MR2 |
| -0.3563329 | 0.6853672 | 0.1031865 | MR2 |
| -0.8453265 | -0.7357813 | -0.4198446 | MR3 |
| -0.0453043 | -0.0588080 | 0.0514924 | MR3 |
| -0.1425694 | -0.1428669 | 0.1159201 | MR3 |
| -0.4840849 | -0.0086084 | -0.7177825 | MR2 |
| 0.0946665 | -0.9830227 | 0.3127275 | MR3 |
| -0.2449772 | -0.1070381 | -0.3615872 | MR2 |
| -0.4055889 | -0.2066661 | -0.6344479 | MR2 |
| 1.9141350 | 1.1447495 | -0.8503421 | MR1 |
| -0.8819515 | -0.7539947 | -0.5299151 | MR3 |
| -0.6088100 | -0.7589018 | -0.5791162 | MR3 |
| -0.3875006 | -0.6174356 | -0.3217773 | MR3 |
| -0.2396999 | -0.6022028 | -0.5971962 | MR1 |
| -0.6937651 | -0.7948259 | 2.4590148 | MR3 |
| 0.5825347 | -0.3969320 | -0.6022342 | MR1 |
| 0.1286559 | -0.3618499 | -0.4197854 | MR1 |
| -0.9057304 | -0.7340345 | -0.4188832 | MR3 |
| -0.6984141 | -0.6291238 | -0.1278735 | MR3 |
| 0.5185966 | 0.9720708 | -0.4334362 | MR2 |
| -0.4783828 | -0.2502371 | -0.4366002 | MR2 |
| 0.6279751 | 1.3743119 | -0.0604470 | MR2 |
| -0.4068809 | -0.9463182 | -0.4246522 | MR1 |
| 1.1170179 | 0.8591049 | -0.6701523 | MR1 |
| 0.5812591 | -0.5721211 | -0.3971179 | MR1 |
| 0.9478023 | -0.7372800 | -0.3699583 | MR1 |
| 0.6130195 | -0.5171601 | 0.0600126 | MR1 |
| 1.3147909 | -0.2059629 | -0.1378673 | MR1 |
| -0.0137440 | 0.1904155 | -0.0136008 | MR2 |
| -0.5892896 | -0.1218455 | -0.2965234 | MR2 |
| 0.4545065 | 0.1776542 | 0.5063414 | MR3 |
| 0.7237581 | -0.3957502 | 0.0910399 | MR1 |
| -0.5699295 | -0.9071261 | -0.1660383 | MR3 |
| -0.5788150 | -0.2426389 | -0.0972306 | MR3 |
| -0.0895472 | -0.2842755 | -0.3541805 | MR1 |
| 1.8829011 | 0.0960242 | -0.0319439 | MR1 |
| -0.3352994 | -0.1798703 | -0.2580813 | MR2 |
| -0.0154998 | -0.3466957 | 0.1796482 | MR3 |
| -0.8172743 | -0.9407225 | -0.1256649 | MR3 |
| 0.8668802 | -0.6254736 | -0.4528465 | MR1 |
| 0.4010778 | -0.7596260 | -0.2304009 | MR1 |
| -0.0714138 | -0.2265309 | -0.3985530 | MR1 |
| 1.0562222 | -0.5087785 | 0.9298282 | MR1 |
| 0.2812909 | -0.2299454 | 0.0760752 | MR1 |
| 1.7480722 | 0.1981629 | 0.5871245 | MR1 |
| 2.4221238 | -0.3465223 | 1.2222252 | MR1 |
| 0.8446953 | 0.2475394 | 0.0520097 | MR1 |
| -0.0951219 | -0.4380578 | -0.1383308 | MR1 |
| 1.0093779 | -0.9543628 | -0.2521802 | MR1 |
| -0.2274832 | -0.6238712 | -0.1637908 | MR3 |
| 0.7782636 | -0.5632695 | 0.6319120 | MR1 |
| 0.1262150 | -0.1895302 | 0.6184334 | MR3 |
| 0.2948640 | 0.1515410 | 1.2712333 | MR3 |
| 2.3103542 | -0.8250749 | 2.1003694 | MR1 |
| -0.1003072 | -0.7572845 | 1.0964998 | MR3 |
| 0.8691671 | -0.1420575 | 1.0413071 | MR3 |
| 0.8404002 | -0.7876054 | 0.5765109 | MR1 |
| 0.9561103 | -0.2276606 | 2.8002004 | MR3 |
| 0.7216891 | -0.3385199 | 1.3366556 | MR3 |
| 0.3970701 | 0.3346541 | 0.5280843 | MR3 |
| 1.0219807 | 0.0694849 | 0.8032418 | MR1 |
| 3.7466863 | -0.5203852 | 1.8487084 | MR1 |
| 2.4720724 | -0.9223499 | 0.7406089 | MR1 |
| 2.1677950 | 0.9227851 | 0.0478785 | MR1 |
| -0.6931873 | -0.2845133 | -0.1162957 | MR3 |
| 1.0707914 | 0.2244252 | 0.2310159 | MR1 |
| -0.1972575 | 0.1457729 | -0.0428793 | MR2 |
| -0.3779183 | -0.5145637 | 0.0309493 | MR3 |
| 0.7697625 | 1.0016966 | 1.5211802 | MR3 |
| -0.4795268 | -0.2658370 | 0.2690625 | MR3 |
| -0.1940969 | -0.5000731 | 0.0859097 | MR3 |
| 0.6709161 | 1.6363516 | 2.1245588 | MR3 |
| 0.4522074 | 0.1777210 | 1.1067403 | MR3 |
| -0.5564833 | -0.9666514 | -0.5073758 | MR3 |
| -0.1274849 | -0.5461634 | 0.0829585 | MR3 |
| -0.5518247 | 0.9082865 | -0.3144294 | MR2 |
| 0.5888204 | 0.3817281 | -0.0636748 | MR1 |
| -0.0966409 | 0.0601297 | 0.1659074 | MR3 |
| -0.6314335 | -0.4335255 | -0.3608555 | MR3 |
| -0.4164280 | -0.2869037 | -0.5322498 | MR2 |
| -1.1907497 | -0.5615732 | -0.4611936 | MR3 |
| -0.1199268 | -0.6545176 | 0.4986514 | MR3 |
| -0.0300377 | 0.4496047 | 0.0514692 | MR2 |
| 0.1203912 | -0.3852759 | 0.4059532 | MR3 |
| -1.5291952 | 6.4688258 | 0.6808807 | MR2 |
| 0.3407989 | 0.3947491 | -0.0014128 | MR2 |
| 0.4934430 | 0.1180693 | 0.9241596 | MR3 |
| -0.6741625 | 0.0858878 | -0.2136856 | MR2 |
| 0.2393436 | -0.2052721 | 0.4144338 | MR3 |
| -0.2723315 | 0.6415275 | 4.9883054 | MR3 |
| -0.8859942 | 1.7716958 | 2.5616129 | MR3 |
| -1.0351382 | -0.8839953 | -0.5161920 | MR3 |
| -0.5975489 | -0.6292137 | -0.1025579 | MR3 |
| -0.2880779 | -0.4809705 | 3.9587162 | MR3 |
| 1.0637374 | 1.2341439 | -0.7029859 | MR2 |
| -0.5546334 | -0.2598207 | -0.5748321 | MR2 |
| -0.4923172 | -0.5174213 | -0.0216993 | MR3 |
| -0.6869757 | 0.0645771 | -0.0005443 | MR2 |
| -1.8045718 | 7.8988815 | 0.1233804 | MR2 |
| -0.2491039 | 0.2102025 | -0.4557714 | MR2 |
| -0.6713179 | -0.4076309 | -0.5522844 | MR2 |
| -0.0095911 | 0.4712679 | -0.3847311 | MR2 |
| -0.0223691 | 0.1733670 | 0.7554614 | MR3 |
| -0.5750973 | 3.4143478 | 0.0205990 | MR2 |
| -0.2603258 | -0.3377156 | -0.5663691 | MR1 |
| -0.2484814 | 0.2224370 | -0.1847181 | MR2 |
| -0.4723640 | -0.3038536 | -0.2581975 | MR3 |
| 0.3614780 | 0.5274274 | -0.1115840 | MR2 |
| -0.5611702 | -0.0524366 | -0.5735648 | MR2 |
| -0.4905256 | -0.0494643 | 0.4030397 | MR3 |
| 0.5398174 | 0.8041059 | -0.6487509 | MR2 |
| -0.2581779 | 0.0451289 | -0.5486750 | MR2 |
| -0.0933303 | 0.3722017 | -0.2281954 | MR2 |
| -0.7167051 | -0.2165924 | -0.4200249 | MR2 |
| -0.7871902 | -0.6426749 | 0.3121102 | MR3 |
| -0.6101046 | -0.2968701 | -0.6367503 | MR2 |
| 0.2121021 | 0.2918378 | 0.5509930 | MR3 |
| -0.7239997 | -0.8430991 | -0.1178101 | MR3 |
| -0.8477597 | -0.0909208 | -0.4824946 | MR2 |
| -0.0300540 | 0.1950799 | -0.0242678 | MR2 |
| -0.6056388 | -0.0927581 | -0.3232018 | MR2 |
| 0.7366014 | 1.3837385 | -0.1821649 | MR2 |
| 1.7386836 | 1.5350040 | -0.5550236 | MR1 |
| -0.9100510 | 0.0800296 | 0.1372773 | MR3 |
| -0.6668953 | -0.5691405 | 0.5284291 | MR3 |
| -0.7107021 | -0.6350926 | 0.5499072 | MR3 |
| -0.2569835 | 1.0775259 | -0.3380717 | MR2 |
| 0.6872389 | 1.6515169 | -0.4594584 | MR2 |
| -0.9638926 | -0.3736237 | -0.2662301 | MR3 |
| -0.6063035 | 0.1191595 | 0.0571289 | MR2 |
| -0.7367727 | -0.1339896 | 0.6710667 | MR3 |
| -0.9579230 | -0.5702865 | 0.9223361 | MR3 |
| -0.7509823 | -0.9743952 | 0.4503113 | MR3 |
| -0.2614695 | 0.0232410 | -0.1945994 | MR2 |
| -0.1138583 | -0.1570644 | 0.4555045 | MR3 |
| 0.3686901 | 0.4783586 | 0.1349106 | MR2 |
| -0.4419724 | -0.2881463 | -0.4939147 | MR2 |
| -0.7479087 | -0.5296224 | -0.4331134 | MR3 |
| 1.1383728 | 1.4997563 | -0.4113166 | MR2 |
| -0.7964850 | -0.6665159 | 0.2229036 | MR3 |
| -0.6942609 | -0.1516090 | -0.3900457 | MR2 |
| -0.5836488 | -0.6199510 | -0.5329672 | MR3 |
| 0.5136656 | 2.0843861 | -0.7368252 | MR2 |
| 2.2336227 | 0.3801604 | 0.3615950 | MR1 |
| -0.6246395 | -0.0758189 | -0.7028974 | MR2 |
| -0.7182380 | -0.5246137 | -0.5962618 | MR2 |
| -0.8003485 | -0.5068291 | -0.4200280 | MR3 |
| -0.6356480 | -0.2842081 | -0.0438185 | MR3 |
| -0.2630800 | -0.0888616 | -0.0681630 | MR3 |
| -0.6386526 | -0.0053009 | -0.5118806 | MR2 |
| 0.6490796 | 0.4944179 | -0.4678989 | MR1 |
| -0.3821492 | -0.7349433 | -0.3852611 | MR1 |
| 0.3913835 | 0.0602158 | -0.6438063 | MR1 |
| 2.1921735 | 0.6816727 | -0.9887817 | MR1 |
| 0.3018942 | 0.6058474 | -0.1769208 | MR2 |
| 2.2061195 | 0.2690747 | -0.5139270 | MR1 |
| 0.0763023 | -0.3397692 | -0.3725075 | MR1 |
| -0.5182625 | 0.0467668 | -0.3550952 | MR2 |
| 0.1700745 | -0.6557965 | -0.6333011 | MR1 |
| 0.5257826 | 0.2164070 | -0.5529576 | MR1 |
| -0.5522442 | -0.6092464 | 0.1887472 | MR3 |
| 0.0379436 | -0.0294007 | 0.0009561 | MR1 |
| -0.5239177 | -1.0435291 | -0.8148104 | MR1 |
| -0.2322537 | -0.4388123 | 0.1680176 | MR3 |
| 3.9462214 | 2.6782029 | 0.2635245 | MR1 |
| -0.3823616 | -1.3623431 | 1.8275739 | MR3 |
| -0.2776017 | -0.8550972 | 0.0088168 | MR3 |
| 0.0015212 | -0.5917970 | -0.5221831 | MR1 |
| -0.1720578 | -0.5467072 | -0.2136313 | MR1 |
| -0.3090690 | -0.5115081 | -0.7569660 | MR1 |
| 1.1512784 | -0.5518301 | 0.5640538 | MR1 |
| 0.5290519 | 0.4737334 | -0.6837971 | MR1 |
| 0.2458078 | 0.4189775 | 0.9215234 | MR3 |
| -0.1996644 | -1.0134231 | 2.1616459 | MR3 |
| 0.0095606 | -0.2772656 | -0.1595063 | MR1 |
| 0.1409053 | -1.1815599 | 0.0227357 | MR1 |
| -0.0384162 | -0.1230969 | -0.2361908 | MR1 |
| 1.2441086 | -0.2347834 | 0.3079196 | MR1 |
| -0.1745638 | -0.5182471 | -0.6356349 | MR1 |
| -0.1097293 | -0.1313849 | 0.0959971 | MR3 |
| -0.4731514 | 0.4109090 | -0.3484272 | MR2 |
| -0.4896922 | -1.1491223 | 2.7968224 | MR3 |
| -0.1485204 | -0.4600702 | 0.0038779 | MR3 |
| -0.1499433 | -0.8032351 | -0.3916093 | MR1 |
| 0.1863467 | -1.0713985 | 0.0329077 | MR1 |
| 1.7929160 | -0.3826485 | -0.5967453 | MR1 |
| 0.8104406 | 2.5263730 | -0.3431188 | MR2 |
| -0.1321465 | 0.1459452 | -0.3282137 | MR2 |
| -0.0001032 | 0.1089340 | -0.6082763 | MR2 |
| -1.1300205 | -0.3258204 | 0.0288609 | MR3 |
| -0.1505416 | -0.2758782 | 0.2045263 | MR3 |
| -0.5251545 | 0.8872251 | 8.8007825 | MR3 |
| -0.8512126 | -0.3886080 | -0.0822090 | MR3 |
| -0.4164706 | 0.1609571 | 7.3514260 | MR3 |
| -0.5759451 | -0.4270808 | -0.2778137 | MR3 |
| -1.1286013 | -0.7539728 | -0.4388792 | MR3 |
| -0.3515784 | 0.1797888 | -0.1373751 | MR2 |
| -0.1433881 | 0.1924909 | 2.2230467 | MR3 |
| -0.6603100 | 0.6909129 | 1.0238736 | MR3 |
| -0.7599532 | 1.8066503 | 0.9163330 | MR2 |
| -0.4340291 | 0.6448290 | 0.6273694 | MR2 |
| -0.6085887 | -0.1388098 | -0.2491734 | MR2 |
| -0.6221689 | -0.3278057 | -0.6958104 | MR2 |
| -0.1958128 | -0.1952106 | 0.1518736 | MR3 |
| -0.1357776 | -0.1512917 | 0.1129693 | MR3 |
| -0.6611733 | -0.4724899 | 0.2293561 | MR3 |
| -1.2246721 | -0.5792883 | 0.3349556 | MR3 |
| -0.4005827 | 0.8415932 | -0.1993459 | MR2 |
| 0.3381614 | 1.5060099 | -0.8926787 | MR2 |
| -1.1162722 | -0.7316743 | -0.4914827 | MR3 |
| -0.4606233 | -0.5643432 | -0.4431495 | MR3 |
| -0.1569905 | 0.4604429 | -0.3249004 | MR2 |
| -0.3863637 | 1.2653365 | -0.2674408 | MR2 |
| 0.2378784 | 0.8505455 | 1.0174546 | MR3 |
| -0.6036365 | -0.5359771 | -0.6168522 | MR2 |
| -0.3887100 | -0.4820988 | 0.1620053 | MR3 |
| -0.8706036 | 0.0765099 | -0.1221859 | MR2 |
| -1.3939527 | 1.2100228 | 0.0447111 | MR2 |
| -0.1977731 | 1.4918446 | -0.5497211 | MR2 |
| -1.1032415 | 0.3667361 | -0.8190350 | MR2 |
| -0.7866845 | 3.4955724 | -0.3813842 | MR2 |
| -1.0857775 | -0.4538314 | -0.0588646 | MR3 |
| 0.0100584 | 1.7499746 | -0.6678818 | MR2 |
| -0.5201353 | -0.0579793 | -0.3966354 | MR2 |
| 0.0496001 | -0.1888492 | -0.6666199 | MR1 |
| -1.3213490 | 0.0702220 | -0.1024905 | MR2 |
| -0.5615563 | 0.7281646 | -0.2378969 | MR2 |
| 0.4158264 | 0.2720363 | -1.0312799 | MR1 |
| -0.8665900 | -0.4042915 | -0.2642658 | MR3 |
| -1.0369352 | 3.1989460 | 0.7973842 | MR2 |
| -0.3407671 | 0.9971294 | -0.9083282 | MR2 |
| -0.9556078 | 2.8078190 | 0.4499527 | MR2 |
| -0.8659092 | 0.1242539 | -0.0562302 | MR2 |
| -0.4135334 | 0.0537588 | -0.8535736 | MR2 |
| -0.8410924 | -0.6358326 | -0.8921011 | MR2 |
| -1.1140889 | -0.1492414 | -0.5899352 | MR2 |
| -0.7040179 | -0.0224744 | -0.9880609 | MR2 |
| -1.2866835 | 0.6977835 | -0.1040257 | MR2 |
| -1.3827310 | 2.7512306 | 1.5018985 | MR2 |
| -1.1420732 | 2.4235500 | 0.7999875 | MR2 |
| -0.8114635 | 0.6087340 | 0.2655096 | MR2 |
| -0.7413819 | 2.3376225 | -0.4279673 | MR2 |
| -0.5830192 | -1.0983433 | -0.4631409 | MR3 |
| -0.9380312 | 1.9160828 | -0.1464056 | MR2 |
| -1.0539815 | -0.1125673 | -0.7825539 | MR2 |
| -0.5969559 | 0.1319212 | 0.2759969 | MR3 |
| -0.5855324 | -0.2874266 | -0.4117792 | MR2 |
| -0.8634694 | -0.0642997 | 0.3079367 | MR3 |
| -0.0892012 | 2.0818549 | -0.4282072 | MR2 |
| -1.2944824 | -0.4943788 | -0.3799764 | MR3 |
| -0.8569657 | -0.2693127 | -0.6126230 | MR2 |
| -0.4419906 | 0.3302242 | 4.0502732 | MR3 |
| -0.1816276 | 0.6549203 | -0.9019410 | MR2 |
| -0.6184108 | -0.6478723 | -0.4955382 | MR3 |
| -0.7690426 | 0.1433578 | -0.5107073 | MR2 |
| 0.5991661 | 0.5607123 | -0.6659591 | MR1 |
| -0.0742518 | 0.3366520 | -0.2374665 | MR2 |
| -1.3972350 | -0.2200077 | 0.0665581 | MR3 |
| -0.9736807 | 1.0943599 | -0.1102189 | MR2 |
| -0.5528531 | 0.4599760 | -0.3228225 | MR2 |
| -0.0974526 | 0.8804731 | -0.3916555 | MR2 |
| -0.6453717 | 0.3126430 | 0.0842114 | MR2 |
| -0.2558817 | -0.6072301 | -0.1679718 | MR3 |
| -0.9197649 | -0.1885489 | -0.6090084 | MR2 |
| 0.2832403 | 1.3221746 | -0.5486807 | MR2 |
| 0.3008081 | 0.8143081 | -0.4476788 | MR2 |
| -0.6815072 | 0.0361392 | -0.3908241 | MR2 |
| 0.6691504 | 3.3649776 | -0.2457160 | MR2 |
| -1.0313253 | 0.2618898 | -0.1056988 | MR2 |
| 0.6986669 | 2.3087502 | -0.4414716 | MR2 |
| -1.1610111 | -1.2241343 | -0.3486476 | MR3 |
| -0.5497612 | -0.9822240 | -0.3442775 | MR3 |
| 0.0890303 | -0.7374852 | -0.8846299 | MR1 |
| -0.5292642 | 0.2894798 | 1.1069114 | MR3 |
| -0.6972670 | -0.7711973 | -0.7827713 | MR1 |
| -0.2586639 | -0.1260834 | -0.4567576 | MR2 |
| 0.4780391 | 0.2912059 | -0.9258342 | MR1 |
| 0.6340053 | -1.0557377 | -0.4311581 | MR1 |
| 0.8195873 | 0.8460390 | 0.9292981 | MR3 |
| 0.8123640 | -0.1162272 | 0.0753595 | MR1 |
| 0.0525879 | 0.1529976 | 0.1086006 | MR2 |
| 0.3894336 | -0.3910573 | -0.3927267 | MR1 |
| -0.9395709 | -1.0041966 | 0.1809076 | MR3 |
| -0.3997114 | 0.3568305 | -0.6345392 | MR2 |
| 0.0210919 | 1.2415189 | -0.8055736 | MR2 |
| 0.1664352 | -0.7872810 | -0.0403167 | MR1 |
| 0.3985694 | 0.0966721 | -0.0117123 | MR1 |
| 0.1269209 | -0.7346039 | -0.4365947 | MR1 |
| -0.5168920 | 0.1798271 | 0.0180594 | MR2 |
| 1.1298070 | 0.2207340 | 0.3354914 | MR1 |
| 0.1180932 | 0.0035832 | -0.8926564 | MR1 |
| 0.3908662 | -0.7744398 | -0.0179717 | MR1 |
| -0.6679958 | -0.4229097 | -0.4154642 | MR3 |
| -0.5682482 | -0.2349497 | 0.0008994 | MR3 |
| 0.1018127 | -0.8019117 | -0.3345540 | MR1 |
| -0.1615954 | -0.9346674 | -0.7259535 | MR1 |
| 0.5245352 | -0.3578485 | -0.3607888 | MR1 |
| -0.8664786 | -0.9180427 | -0.1242794 | MR3 |
| 0.1796114 | -0.8412609 | 0.1062423 | MR1 |
| 0.5582890 | 0.3309434 | 0.2265501 | MR1 |
| 0.7504296 | 0.3216100 | -0.8160629 | MR1 |
| -0.8445016 | -0.0386561 | -0.0627013 | MR2 |
| -0.3787751 | -0.2753616 | 1.2250517 | MR3 |
| 0.1671814 | -0.1104891 | 0.1210681 | MR1 |
| -0.9496812 | 0.0640661 | -0.8410155 | MR2 |
| -0.4962602 | -0.4639897 | -0.7977540 | MR2 |
| -0.7143047 | 0.0951926 | 1.5860021 | MR3 |
| -0.7533885 | -0.9401264 | -0.4796811 | MR3 |
| -0.7166842 | -0.8646194 | 0.1013292 | MR3 |
| 0.2162779 | -0.4746694 | -0.7432182 | MR1 |
| -0.6182079 | -0.3890225 | 0.0007661 | MR3 |
| -0.8729683 | -0.1758211 | 0.6448502 | MR3 |
| -0.5808042 | -0.4126343 | -0.6038459 | MR2 |
| -0.5143599 | -0.1896185 | -0.3865641 | MR2 |
| -1.1345970 | -0.6993779 | -0.2478149 | MR3 |
| -0.2810240 | -0.2674838 | 0.5819758 | MR3 |
| -0.9493428 | -0.5698946 | -0.5016026 | MR3 |
| 0.0274796 | 0.9025948 | 1.4218189 | MR3 |
| -1.0107792 | 0.4255178 | -0.3945677 | MR2 |
| -0.2901718 | 0.9898824 | -0.5287472 | MR2 |
| -0.8917733 | -1.1150342 | -0.3569310 | MR3 |
| -0.4122739 | -0.2056138 | -0.1058017 | MR3 |
| -0.5134281 | -0.7691100 | 0.1522760 | MR3 |
| 1.1876186 | 1.0612328 | 0.0985591 | MR1 |
| -0.4570197 | -0.2319029 | 4.0054156 | MR3 |
| -1.2570058 | 0.0604841 | 0.1136132 | MR3 |
| -1.2755569 | -0.6657473 | -0.4020904 | MR3 |
| -0.2258289 | -0.3512983 | -0.6181290 | MR1 |
| -1.1605933 | -0.2581708 | -0.0922461 | MR3 |
| -0.5993678 | 0.3144513 | -0.8256399 | MR2 |
| 0.5985531 | 2.1017443 | 0.0585061 | MR2 |
| -0.9413155 | 0.0691136 | -0.4291663 | MR2 |
| -1.0924473 | -0.6315343 | -0.8828704 | MR2 |
| -1.0034272 | -0.8213221 | -0.2940829 | MR3 |
| -0.9354530 | -0.4760207 | -0.9006219 | MR2 |
| -1.0973423 | 0.7308503 | 0.2363849 | MR2 |
| -0.9216153 | -0.5216837 | -0.3994497 | MR3 |
| -0.3540298 | 1.0688919 | -0.5953494 | MR2 |
| 0.4552506 | 1.5054561 | -0.5385139 | MR2 |
| -1.1257309 | 0.0925854 | 0.3396387 | MR3 |
| -1.0919093 | 0.4385398 | 5.3103265 | MR3 |
| -0.3316389 | 1.0247186 | -0.4805893 | MR2 |
| -0.8712678 | -0.0476240 | -0.6820480 | MR2 |
| -0.1022288 | -0.2434787 | -0.5473678 | MR1 |
| -0.2333178 | 0.4707045 | 0.4141105 | MR2 |
| 0.2236249 | 1.1382155 | -0.9563246 | MR2 |
| -0.8311771 | -0.5563502 | -0.0654780 | MR3 |
| -0.8992964 | -0.1447900 | -0.4487358 | MR2 |
| -0.2256448 | 0.5897143 | -0.5025268 | MR2 |
| -0.2536727 | 0.5241444 | 0.1969867 | MR2 |
| -1.2452084 | -0.6140507 | -0.3521655 | MR3 |
| -0.5251419 | 0.2478437 | -0.7249974 | MR2 |
| -0.7245745 | -0.9185332 | -0.3771058 | MR3 |
| -0.7112391 | -0.0608165 | -0.6559277 | MR2 |
| -0.9901076 | 0.1820206 | 0.2199504 | MR3 |
| -0.0526492 | 0.7126386 | -0.2607844 | MR2 |
| -1.1871033 | -0.2808652 | -0.1974240 | MR3 |
| -1.1052735 | -0.0692250 | -0.9235143 | MR2 |
| -0.7417229 | -0.2410923 | -0.7539784 | MR2 |
| -0.5828891 | -0.4573417 | -0.4382042 | MR3 |
| -0.7942265 | -0.0208897 | 0.1509365 | MR3 |
| -0.8639451 | -0.1913489 | -0.2115774 | MR2 |
| -0.5481511 | -0.2182146 | -0.8668657 | MR2 |
| -0.9473471 | -0.0984660 | 0.0142281 | MR3 |
| -1.0723975 | -0.4041846 | -0.3576521 | MR3 |
| -1.0992925 | -0.2733968 | -0.3189757 | MR2 |
| -0.6230211 | 0.0178861 | -0.3716321 | MR2 |
| -0.3711339 | 0.0281853 | -0.9751057 | MR2 |
| -1.2039141 | 0.1184809 | -0.2203242 | MR2 |
| 0.0419752 | 1.4971126 | 0.1881315 | MR2 |
| 0.9047104 | -0.9010946 | -0.3912974 | MR1 |
| -2.3897140 | 8.7889910 | 1.6337526 | MR2 |
| -1.1394055 | 0.1137202 | -0.1952212 | MR2 |
| -0.3239849 | 0.8572591 | -0.8207892 | MR2 |
| -1.0446108 | -0.3172831 | -0.2529293 | MR3 |
| -0.8214552 | 0.4449314 | -0.0825795 | MR2 |
| -0.5531485 | 0.1105514 | -0.0571077 | MR2 |
| -0.7327056 | 0.6814485 | 1.3898716 | MR3 |
| -0.4853038 | 0.1265499 | -1.0109955 | MR2 |
| -0.2044355 | -1.0982881 | 2.2877495 | MR3 |
| -1.0079662 | -0.2692160 | -0.1185834 | MR3 |
| -0.5553415 | 1.2681913 | 0.1173336 | MR2 |
| -0.6406597 | 1.8602343 | -0.1360138 | MR2 |
| -0.5747073 | 1.1511972 | -0.6678180 | MR2 |
| -0.2957165 | -0.5564438 | 1.3586057 | MR3 |
| -0.3618281 | -0.2294092 | -0.2577325 | MR2 |
| -0.4571332 | -0.1309541 | -1.0442294 | MR2 |
| -0.5325381 | -0.5438308 | -0.0126916 | MR3 |
| -0.8139444 | 0.1196772 | 0.7299056 | MR3 |
| -0.2182981 | 0.3163416 | -0.3163155 | MR2 |
| -0.8788343 | -0.4582744 | 0.0232076 | MR3 |
| -0.3874267 | 0.2847144 | 0.1118690 | MR2 |
| 0.1015213 | -0.6122387 | -0.5162445 | MR1 |
| 0.0829123 | -0.4067746 | -0.0432755 | MR1 |
| -0.7127719 | -0.0634838 | -0.4587994 | MR2 |
| -0.7840505 | -0.3205866 | -0.3957931 | MR2 |
| 0.3159858 | -0.7476342 | -0.3241250 | MR1 |
| -0.6022923 | -0.1431967 | -0.8167778 | MR2 |
| 0.3619877 | -1.1625312 | -0.9011732 | MR1 |
| 0.1675582 | -0.5662077 | -0.5087182 | MR1 |
| 0.0540215 | -0.5403880 | 0.1443664 | MR3 |
| 1.3680908 | 0.0118387 | -0.7978029 | MR1 |
| -0.7036948 | -0.0444695 | 0.6156921 | MR3 |
| -0.6568407 | -0.7345504 | -0.1314824 | MR3 |
| 0.0509130 | -0.6438263 | -0.6151874 | MR1 |
| 1.0886326 | -0.8233934 | -0.7184992 | MR1 |
| 0.8900928 | -0.6504151 | -0.4196888 | MR1 |
| -0.0018729 | -1.0211490 | 0.3462723 | MR3 |
| 1.2255653 | -1.5051635 | -0.1657137 | MR1 |
| 0.7140234 | -0.3758688 | -0.4243409 | MR1 |
| 0.0479973 | -0.9421957 | -0.4151628 | MR1 |
| 0.4856884 | -0.1752931 | 0.0571939 | MR1 |
| -0.0579543 | 0.3431011 | -0.8935474 | MR2 |
| -0.6563109 | -0.4230745 | -0.7161283 | MR2 |
| 0.0992763 | -0.1931487 | -0.3318996 | MR1 |
| -0.3978011 | -0.6938737 | -0.6172264 | MR1 |
| -0.0661869 | 0.2438537 | 0.2607997 | MR3 |
| -0.1577048 | -0.4796741 | 0.0413021 | MR3 |
| 0.4230178 | -0.7782144 | -0.0889169 | MR1 |
| 1.4000645 | 1.1474869 | -0.4249125 | MR1 |
| -0.7581524 | -0.3885536 | -0.5252251 | MR2 |
| -0.6250532 | 0.2711556 | -0.0158182 | MR2 |
| -0.4529672 | -0.3039790 | 0.5699610 | MR3 |
| -0.8755071 | -1.0827340 | -0.1830954 | MR3 |
| -0.3053638 | -0.9350218 | -0.5822589 | MR1 |
| -0.0433133 | -0.0903196 | -0.3279478 | MR1 |
| -0.2813722 | -0.8838322 | -0.5294480 | MR1 |
| -0.7984354 | -0.5676808 | -0.2456271 | MR3 |
| -0.4193455 | -0.1910420 | 0.1159450 | MR3 |
| -1.0700681 | -0.4582830 | -0.2947169 | MR3 |
| -0.1400873 | -0.1118479 | -0.4728551 | MR2 |
| -1.1223347 | 0.1451876 | 0.2753565 | MR3 |
| -0.3499401 | -0.8220289 | 0.8909544 | MR3 |
| -0.4974500 | 0.1835915 | 0.1650293 | MR2 |
| -0.5121940 | -0.3680601 | 0.5311614 | MR3 |
| -0.2250349 | -0.0192928 | 0.3480621 | MR3 |
| -0.5147260 | -1.2743509 | 0.0801025 | MR3 |
| -0.4789884 | -0.5560672 | -0.5409548 | MR1 |
| -0.6600478 | -0.0896841 | 0.2761180 | MR3 |
| -0.0982587 | -0.3029336 | 0.2385141 | MR3 |
| -0.5094703 | 3.0512379 | 5.7563869 | MR3 |
| -0.7664082 | -0.2060160 | -0.3457530 | MR2 |
| -0.6972460 | -0.9219253 | -0.4402427 | MR3 |
| 0.7084541 | 0.8147353 | -0.7015627 | MR2 |
| -0.2511796 | -0.6325017 | -0.8467577 | MR1 |
| -0.8149299 | -0.6336528 | -0.6191591 | MR3 |
| -0.5741802 | 3.2578766 | 0.3814191 | MR2 |
| -0.9595433 | 0.9041255 | 0.1752967 | MR2 |
| -0.6102317 | 0.2886808 | -0.3547575 | MR2 |
| -0.6749851 | 2.1695569 | -1.0168479 | MR2 |
| 0.6848148 | -1.1938935 | -0.6024127 | MR1 |
| 1.0920726 | -0.5788889 | -0.7221029 | MR1 |
| -0.2051813 | -1.0537039 | -1.2496361 | MR1 |
| 0.2145003 | 0.2465911 | -1.2699997 | MR2 |
| 0.1770566 | 1.1959258 | 0.4294363 | MR2 |
| 0.0069436 | 1.0263651 | -0.2360369 | MR2 |
| 0.4262603 | -0.0931654 | -1.1311904 | MR1 |
| -0.6561684 | -0.9831610 | -0.8024747 | MR1 |
| -0.8143609 | -0.7028815 | -0.0821071 | MR3 |
| -0.5737506 | 0.3209525 | 0.0122874 | MR2 |
| -0.7147055 | 0.1765779 | 0.5420513 | MR3 |
| -0.8253046 | 0.9183296 | -0.7962719 | MR2 |
| -1.2350067 | 1.5439055 | -0.6205585 | MR2 |
| -0.1379146 | 1.3676130 | -1.0244250 | MR2 |
| -0.3792393 | -0.0090357 | -0.9008670 | MR2 |
| -0.7708184 | 0.4908518 | 0.5907342 | MR3 |
| 0.8959425 | 1.3772582 | -0.9159293 | MR2 |
| -0.1454651 | 1.3852921 | 0.9799499 | MR2 |
| -1.3321508 | -0.2006044 | -0.0085502 | MR3 |
| -0.0409552 | 1.6214165 | -0.8435577 | MR2 |
| -0.7270276 | 0.1920919 | -0.5606008 | MR2 |
| -0.1174810 | 1.7055662 | 0.3335779 | MR2 |
| -1.2498492 | 0.2211910 | 0.0046501 | MR2 |
| -1.0970094 | -0.1273818 | 6.1493362 | MR3 |
| -0.7108537 | -0.4924562 | 3.2306056 | MR3 |
| -0.1049207 | 0.4871469 | -0.6765685 | MR2 |
| -0.3150224 | -0.2959662 | -0.6208410 | MR2 |
| -0.5407781 | -0.6870277 | 0.4108801 | MR3 |
| 1.4689092 | 1.2233878 | -0.4768299 | MR1 |
| -0.3108506 | 0.3559553 | -0.9251810 | MR2 |
| -0.2418798 | -0.2132813 | -0.3664917 | MR2 |
| -0.2212713 | 0.5996921 | -0.6295044 | MR2 |
| -1.7338532 | 4.3363535 | 0.0015901 | MR2 |
| -1.0813262 | -0.4731526 | 0.4218633 | MR3 |
| -0.5930276 | 1.4299654 | -0.4893355 | MR2 |
| 0.7806854 | 0.9734987 | -0.6978176 | MR2 |
| -0.4243243 | 0.5839121 | -0.8620330 | MR2 |
| -1.0010024 | -0.0075732 | 0.8172457 | MR3 |
| -0.2603553 | -0.0109133 | -0.7110095 | MR2 |
| -1.0158508 | -0.3241436 | -0.5683718 | MR2 |
| 0.0297618 | -0.0040565 | 0.0028653 | MR1 |
| -0.9465294 | -0.5559364 | -0.4781119 | MR3 |
| -0.8456300 | 0.0393195 | 0.5849587 | MR3 |
| -1.1927697 | -0.6096475 | 0.0691575 | MR3 |
| -1.0868169 | -0.2008635 | -0.1450368 | MR3 |
| -0.0936877 | -0.1069702 | 0.0842286 | MR3 |
| -0.2562464 | 0.4973167 | -0.6797866 | MR2 |
| -0.0935783 | 0.0888178 | -0.3658006 | MR2 |
| -1.3279436 | 0.0474372 | 0.2312443 | MR3 |
| -0.4106439 | -0.6926586 | 0.3223236 | MR3 |
| -0.0861263 | -0.5230949 | 1.6459884 | MR3 |
| -0.3305412 | -0.0825310 | -0.4930620 | MR2 |
| -0.3512888 | 0.4941761 | -0.5788111 | MR2 |
| -0.4358465 | 0.1205028 | -0.8588698 | MR2 |
| -0.5670258 | 0.6189890 | 0.1173626 | MR2 |
| -0.5519591 | -0.3474835 | -0.7512300 | MR2 |
| -0.4335575 | -0.3979723 | -0.2308991 | MR3 |
| -0.1219039 | 1.4361334 | 1.0477080 | MR2 |
| -0.6041916 | -0.0636204 | -0.3363625 | MR2 |
| 0.0589340 | 0.0137419 | -0.0156554 | MR1 |
| -1.0227108 | -1.0872137 | -0.1353724 | MR3 |
| 0.2044725 | -0.5160457 | -0.7285424 | MR1 |
| -0.9599369 | -0.5566707 | -0.0053378 | MR3 |
| 2.1273341 | 0.2360225 | -0.5569595 | MR1 |
| -0.7788776 | -1.1226875 | -0.1339539 | MR3 |
| -0.8204992 | 0.3338060 | -0.5880195 | MR2 |
| -0.5737594 | -0.0247095 | 0.1133780 | MR3 |
| 0.4231114 | -0.5851061 | -0.4754964 | MR1 |
| -0.6058951 | -0.9711655 | 0.0916745 | MR3 |
| -0.5848098 | 0.0778334 | -0.2520311 | MR2 |
| 0.1376571 | -0.5062786 | -0.4926904 | MR1 |
| -0.3850137 | -0.5418023 | -0.1271780 | MR3 |
| -1.0230841 | -0.8200198 | 0.5379026 | MR3 |
| 1.1468059 | -0.2621025 | -0.0522270 | MR1 |
| -0.6235918 | -1.6834906 | 2.4791385 | MR3 |
| -0.8198157 | -0.9669911 | -0.2964077 | MR3 |
| -0.3580613 | -0.9139771 | 1.1836230 | MR3 |
| -0.0972338 | 0.3252576 | -0.0918687 | MR2 |
| -0.1627650 | -0.1460506 | -0.5439976 | MR2 |
| 0.2622780 | -1.3104409 | -0.2755282 | MR1 |
| -0.1099317 | -0.6884313 | -0.4421060 | MR1 |
| -0.0381546 | -0.3408859 | 0.1823348 | MR3 |
| 0.7123492 | 0.3398816 | -0.3382229 | MR1 |
| 1.5250097 | -0.8305444 | -0.3782219 | MR1 |
| -0.0185770 | -0.9853536 | -0.4044916 | MR1 |
| -0.8840227 | -0.9641181 | 1.2352069 | MR3 |
| -0.0864887 | 0.6136389 | 1.0526131 | MR3 |
| 0.7723600 | -0.3243385 | 1.4948480 | MR3 |
| 2.0252272 | -0.6761340 | 0.1928634 | MR1 |
| 2.3035453 | 1.4315972 | -0.8268883 | MR1 |
| -0.0105703 | -0.4487548 | -0.4301255 | MR1 |
| 0.7666346 | -0.7078051 | -0.4433114 | MR1 |
| 0.4987887 | -0.6652850 | -0.3136437 | MR1 |
| -0.7706531 | -0.4586352 | -0.6715432 | MR2 |
| -0.2644888 | -1.1076953 | -0.3159514 | MR1 |
| -0.4635745 | -0.4689144 | -0.0886921 | MR3 |
| -0.5108487 | -0.3908129 | -0.4161846 | MR2 |
| -0.4125132 | -0.8880089 | -0.5716112 | MR1 |
| -0.2707010 | -0.8345025 | -0.0107304 | MR3 |
| 0.9001773 | -0.9882452 | 0.2240771 | MR1 |
| 0.3073828 | -1.4452187 | 0.4603304 | MR3 |
| 0.0235657 | 0.3962700 | 0.2190011 | MR2 |
| -0.7690786 | -1.1289550 | 0.0754598 | MR3 |
| 0.5636260 | -0.4449240 | -0.0228834 | MR1 |
| -0.5202355 | -0.8171045 | 0.1645487 | MR3 |
| 0.0803989 | -0.0633024 | 0.8753476 | MR3 |
| -0.7584444 | -0.7554195 | 2.2938990 | MR3 |
| -0.6048950 | -0.8958718 | 2.6304418 | MR3 |
| -0.0355754 | 0.0831579 | -0.2660020 | MR2 |
| -0.0907380 | -0.0925737 | -0.7775577 | MR1 |
| -0.0190994 | -0.8218593 | 0.5105178 | MR3 |
| -0.8490334 | -0.2243330 | 0.6478017 | MR3 |
| -1.0829263 | -0.4469538 | -0.0453849 | MR3 |
| -1.1280980 | -0.2776348 | 3.3375469 | MR3 |
| -0.8401364 | -0.7258659 | 0.4073790 | MR3 |
| -0.9034115 | 0.7457321 | -0.7044293 | MR2 |
| -0.6674330 | -0.5402265 | -0.6615114 | MR2 |
| -0.4670659 | -0.2916108 | -0.3446936 | MR2 |
| -0.3967094 | -1.0492825 | 0.1313441 | MR3 |
| -1.0844562 | -0.0968413 | -0.4690062 | MR2 |
| -0.7888606 | -0.7129837 | 0.5853575 | MR3 |
| -0.9041157 | -1.2680586 | -0.3629005 | MR3 |
| -0.5135259 | 0.1917535 | 0.6622811 | MR3 |
| -0.7588962 | 0.5086753 | -0.1618549 | MR2 |
| -0.9851330 | 1.2058653 | 0.0599963 | MR2 |
| -0.5477924 | 0.3297233 | -0.3907021 | MR2 |
| -0.7360724 | -0.5270074 | -0.4282366 | MR3 |
| -0.3775558 | -0.4275496 | -0.6528845 | MR1 |
| -0.0080696 | 0.1614023 | -0.6862866 | MR2 |
| -0.6842766 | 0.9254762 | 1.7788220 | MR3 |
| -0.1883554 | -0.1371336 | -0.4066902 | MR2 |
| -0.0456297 | -0.0694677 | 0.0528318 | MR3 |
| -0.4999201 | -0.2944899 | 0.0079380 | MR3 |
| -0.8219232 | -0.0935319 | -0.1039839 | MR2 |
| -0.2716116 | -0.5036760 | -0.5670476 | MR1 |
| -0.7345844 | -0.6539050 | -0.7626773 | MR2 |
| 0.0718023 | 0.0643843 | -0.0284591 | MR1 |
| 0.2286859 | 0.2296023 | -0.1355871 | MR2 |
| -0.2084814 | 1.6377383 | -0.6815411 | MR2 |
| -0.0905572 | 0.0728217 | 0.0587635 | MR2 |
| -0.4008610 | -0.9246766 | -0.7601572 | MR1 |
| -0.2399118 | 0.5661512 | 0.0312686 | MR2 |
| 0.0019819 | -0.3508475 | -0.5534886 | MR1 |
| -0.1581737 | 1.5225474 | -0.2338131 | MR2 |
| -0.3771097 | 3.2830089 | -0.4344141 | MR2 |
| -0.0787728 | -0.0976704 | 0.0747379 | MR3 |
| -0.1749924 | 0.5696376 | -0.9614572 | MR2 |
| 0.4830693 | 1.1130607 | 0.2753721 | MR2 |
| -0.4480813 | 0.4214026 | -0.5822584 | MR2 |
| -0.0911536 | -0.4878375 | -0.3195304 | MR1 |
| -0.3971517 | -0.7167260 | -0.1221555 | MR3 |
| -0.9749050 | -0.0579235 | -0.4300923 | MR2 |
| -0.1294638 | 0.2164683 | -0.7073834 | MR2 |
| -0.5232110 | -0.3144365 | 0.9858810 | MR3 |
| -0.3602378 | 0.1908459 | -0.1731460 | MR2 |
| -0.3123762 | 1.7554728 | -0.3035835 | MR2 |
| -0.6617725 | -1.1942990 | -0.1881562 | MR3 |
| -0.4100746 | -0.9220246 | -0.3330486 | MR3 |
| -0.7698442 | 0.3160878 | -0.0493924 | MR2 |
| -0.8661284 | -0.1633838 | -0.5399745 | MR2 |
| -1.2185213 | -0.3044585 | -0.1064239 | MR3 |
| -0.5443314 | -0.6313586 | -0.4766675 | MR3 |
| -0.2143027 | 0.3711279 | -0.2735811 | MR2 |
| 0.0722588 | -0.6888282 | 0.2276015 | MR3 |
| -0.6653328 | -0.7233718 | -0.2800998 | MR3 |
| -0.3302988 | 2.8406558 | 0.4008092 | MR2 |
| -0.7617739 | -0.7482802 | -0.3347320 | MR3 |
| -0.7236502 | 0.2107789 | -0.6916028 | MR2 |
| -0.8271158 | 0.6051999 | -0.4951979 | MR2 |
| -1.2631856 | -0.8872937 | -0.2874304 | MR3 |
| -1.4056786 | -0.3828130 | 3.6439085 | MR3 |
| -0.6561763 | -1.1919726 | -0.5966421 | MR3 |
| -0.5343242 | -0.6825320 | -0.2365506 | MR3 |
| -0.2646712 | -0.4183841 | -0.5990272 | MR1 |
| -0.6490295 | -1.1728693 | -0.3683732 | MR3 |
| -0.3007383 | 0.0412676 | -1.0318305 | MR2 |
| 0.0172150 | 0.2981445 | -0.5609501 | MR2 |
| -1.3394499 | -0.4196013 | 0.1321899 | MR3 |
| -1.0985654 | 0.0522089 | -0.2719795 | MR2 |
| -0.6110674 | 0.4771421 | 0.1411815 | MR2 |
| -0.7346754 | -0.3472811 | -0.2313898 | MR3 |
| -0.5658373 | -0.9626572 | -0.3720028 | MR3 |
| 0.8446927 | 0.6980086 | -0.5924419 | MR1 |
| -0.5651782 | 0.1586013 | 0.6004054 | MR3 |
| -1.0475880 | -0.8223389 | 4.3934250 | MR3 |
| -0.2681592 | 1.1205868 | 0.7777050 | MR2 |
| -0.6146311 | -0.4371207 | -0.5831938 | MR2 |
| 0.3548021 | -0.1781395 | -0.6841926 | MR1 |
| -1.1548031 | -0.0643735 | -0.2019030 | MR2 |
| -1.0846475 | -0.4547449 | -0.1647718 | MR3 |
| -0.6257415 | -0.0069499 | -0.4827148 | MR2 |
| -1.0354958 | -0.4283108 | 0.1877239 | MR3 |
| -0.4722060 | -1.1159107 | -0.0085501 | MR3 |
| -1.0667982 | 0.6774868 | -0.1750077 | MR2 |
| -0.8551728 | -0.6020359 | 0.8162270 | MR3 |
| -0.1010936 | 0.0153990 | 0.1171803 | MR3 |
| 2.8042590 | 0.2188109 | -0.5577446 | MR1 |
| 0.6443465 | -0.5057973 | 0.1362443 | MR1 |
| 0.1798766 | 0.2473856 | 0.3684131 | MR3 |
| -0.7506269 | -0.6455706 | -0.7476164 | MR2 |
| 0.0074180 | -0.1884924 | -0.6259377 | MR1 |
| 0.9081200 | -0.5247182 | -0.3140219 | MR1 |
| 0.4058742 | 1.3114118 | 1.6443745 | MR3 |
| -0.7869525 | -1.0090170 | 1.0238601 | MR3 |
| 0.3517481 | -1.0941926 | -0.3642128 | MR1 |
| -0.5130203 | -1.1489682 | -0.4768475 | MR3 |
| -0.7720570 | -1.2122234 | 0.3636060 | MR3 |
| -0.4279593 | -0.8557941 | -0.5163770 | MR1 |
| 0.0530853 | -0.8537101 | -0.2810657 | MR1 |
| 0.2459291 | 0.4040516 | -0.5129087 | MR2 |
| -0.8055914 | -0.5512830 | 0.2493422 | MR3 |
| -0.1029910 | -1.0924682 | -0.4882713 | MR1 |
| 0.6660240 | -0.7313207 | -0.2594983 | MR1 |
| -0.2796650 | -0.7721497 | 0.1297033 | MR3 |
| -0.1358032 | -0.1683227 | 1.0263921 | MR3 |
| -0.7044104 | -0.1960287 | -0.4961629 | MR2 |
| 0.7369550 | 0.7283248 | -0.5935518 | MR1 |
| -0.7809892 | -0.5189501 | -0.4114421 | MR3 |
| -0.2846904 | -0.9774096 | -0.4477862 | MR1 |
| -0.9643944 | -0.2373787 | -0.3404013 | MR2 |
| -0.9929569 | 0.1919658 | -0.4084467 | MR2 |
| -0.4147958 | 0.5642268 | 0.4626391 | MR2 |
| -1.0643713 | 0.0055112 | 0.1899322 | MR3 |
| -1.1727989 | -0.2499761 | 0.6133901 | MR3 |
| -0.7179500 | -0.4386473 | -0.6348845 | MR2 |
| -0.8958939 | -1.2861646 | -0.7482955 | MR3 |
| -0.7322747 | -1.7982377 | 13.6976009 | MR3 |
| -0.3163309 | -0.3681656 | -0.4651614 | MR1 |
| -0.9225809 | -0.2528515 | 0.0549442 | MR3 |
| -1.0611162 | -1.1243773 | -0.5178523 | MR3 |
| -0.2749664 | -0.9473949 | -0.4111990 | MR1 |
| 0.5087940 | 1.3101170 | -0.3124190 | MR2 |
| -0.7108336 | -1.2052149 | -0.1573144 | MR3 |
| -0.3950552 | -0.9362000 | 3.5670734 | MR3 |
| -0.6498714 | -0.8294318 | 2.0060665 | MR3 |
| -0.8218672 | -0.7856838 | -0.1790562 | MR3 |
| -0.5669438 | -0.4531847 | -0.2130888 | MR3 |
| -0.9337908 | 1.6302873 | -0.4978775 | MR2 |
| -0.4732510 | 0.1414143 | -0.1896123 | MR2 |
| -0.7293355 | -0.5946153 | -0.4077981 | MR3 |
| -0.1368324 | -0.1549124 | 5.9531284 | MR3 |
| -0.9592217 | -0.3778684 | -0.3196026 | MR3 |
| -0.6490755 | -0.7036140 | -0.4262930 | MR3 |
| -0.1461405 | -0.0175107 | 0.0717356 | MR3 |
| -0.1677209 | -0.0699731 | -0.1751359 | MR2 |
| -0.7284760 | 0.7534271 | -0.6892293 | MR2 |
| -0.4036824 | -0.1524026 | 0.1131148 | MR3 |
| -0.4615044 | 1.0349226 | -0.6294180 | MR2 |
| -0.3124259 | 1.2687829 | 0.0711831 | MR2 |
| -0.7411685 | -0.2975270 | 0.1613333 | MR3 |
| -0.5499302 | 0.3969163 | 0.2292556 | MR2 |
| -0.5422027 | -0.2388488 | -0.9914562 | MR2 |
| -0.2759000 | -0.5770390 | -0.3348578 | MR1 |
| -0.3961931 | 1.4553775 | -0.1887066 | MR2 |
| -0.0766055 | -0.0227336 | -0.2512505 | MR2 |
| -0.6873514 | 0.1727206 | 0.2466783 | MR3 |
| -0.8733115 | 1.2563800 | 0.5861007 | MR2 |
| -0.0317995 | 0.0041165 | -0.8506112 | MR2 |
| -1.0741870 | -0.9695124 | -0.5426867 | MR3 |
| -0.2224276 | -0.9183435 | -0.4228661 | MR1 |
| -0.1394345 | -0.0139766 | -0.4926637 | MR2 |
| -0.3980606 | 0.0869997 | 0.2699113 | MR3 |
| -0.7821795 | -1.1842021 | -0.3200820 | MR3 |
| -0.4322197 | 0.3777025 | -0.2775068 | MR2 |
| -0.1886159 | 1.0862288 | -0.7310323 | MR2 |
| -0.1239381 | -0.2639592 | -0.8442132 | MR1 |
| -0.3146950 | -0.5895244 | -0.7233452 | MR1 |
| -0.9941256 | -0.1472027 | -0.1254009 | MR3 |
| -0.8830363 | -0.4376730 | -0.7521494 | MR2 |
| -0.5617927 | -0.3537958 | -0.3339589 | MR3 |
| -0.4947488 | -0.7963660 | -0.9300862 | MR1 |
| -1.1706590 | -0.5014468 | 0.9132331 | MR3 |
| -1.2975719 | -0.3481347 | -0.5192587 | MR2 |
| -0.9767753 | -0.4370971 | -0.9683655 | MR2 |
| -1.0959860 | 0.4506416 | -0.1324284 | MR2 |
| 0.6143130 | 0.8092216 | -0.7731869 | MR2 |
| -0.2681979 | -0.8139357 | 0.5167577 | MR3 |
| -0.6966731 | -0.8180578 | -0.5577791 | MR3 |
| -1.1772341 | 0.2586295 | 0.3516751 | MR3 |
| -0.4706045 | -0.0363026 | -0.4557712 | MR2 |
| -0.6976780 | 0.3927212 | -0.3000132 | MR2 |
| -0.9449388 | -1.1400894 | 0.2123225 | MR3 |
| -0.7059354 | -0.3630221 | 0.5345395 | MR3 |
| -0.5277002 | 0.7109439 | -0.8129243 | MR2 |
| -1.1379215 | -0.7422654 | -0.1125592 | MR3 |
| -0.9357997 | -0.2744292 | -0.6986634 | MR2 |
| -0.8293432 | -0.1327143 | -0.6455610 | MR2 |
| -1.0394213 | -0.6434912 | -0.4575731 | MR3 |
| -1.1367833 | -0.5695075 | -0.2531394 | MR3 |
| -0.2836723 | 0.4900244 | -0.5946472 | MR2 |
| 1.7044828 | 2.2395612 | -0.8822518 | MR2 |
| -0.5737762 | -0.9750622 | 0.8215575 | MR3 |
| -0.0031482 | -0.4323106 | 0.4836113 | MR3 |
| -0.5526019 | -0.2619633 | 0.2027503 | MR3 |
| -0.7462665 | -0.7596980 | 0.6287715 | MR3 |
| -0.8075750 | -0.1058168 | 1.2931077 | MR3 |
| -0.2924270 | -0.0114802 | 0.2133922 | MR3 |
| -0.5469644 | -0.2509861 | -0.6123954 | MR2 |
| -0.5463751 | -0.5576759 | -0.4253122 | MR3 |
| -0.4984889 | -0.3218890 | 0.0353720 | MR3 |
| 1.1971087 | 0.8965839 | -0.4003693 | MR1 |
| -0.6174029 | -0.4492717 | -0.6164126 | MR2 |
| -0.8743447 | -0.0838477 | -0.5366796 | MR2 |
| 0.6570066 | -0.1263193 | -0.1554178 | MR1 |
| -0.5506899 | -0.2837092 | -0.6856630 | MR2 |
| 0.0997294 | -1.0871860 | -0.7073075 | MR1 |
| -0.9064675 | -0.5612589 | -0.1632195 | MR3 |
| -0.8054708 | -0.1320078 | 0.1263702 | MR3 |
| -0.7252666 | -0.5142194 | -0.4850341 | MR3 |
| -0.8855629 | -0.1302667 | -0.3017862 | MR2 |
| -0.7029732 | -0.2042943 | -0.3501229 | MR2 |
| 2.0713605 | 0.0401420 | 0.3479842 | MR1 |
| -0.9107427 | -0.7384042 | -0.2664497 | MR3 |
| -0.1969714 | -0.5281011 | -0.6516813 | MR1 |
| -0.0855104 | 0.1062078 | -0.6147346 | MR2 |
| -0.5155914 | 0.3582580 | -0.6458996 | MR2 |
| 0.0934059 | -0.1350080 | -0.5984857 | MR1 |
| -0.3666219 | -0.8193410 | -0.1353239 | MR3 |
| -0.5411815 | -0.3781702 | -0.4361143 | MR2 |
| 0.1946964 | 0.2185054 | -0.8113513 | MR2 |
| 0.6139677 | -0.5815633 | -0.2314331 | MR1 |
| -0.3087861 | -0.3388645 | -0.4207362 | MR1 |
| -0.3966363 | -0.8710518 | -0.5983809 | MR1 |
| -1.1062406 | -0.4489055 | -0.1614480 | MR3 |
| 0.1430667 | -0.2599751 | -0.4400106 | MR1 |
| 0.1843410 | -0.1097954 | 0.1176128 | MR1 |
| -0.2932372 | -1.0172667 | -0.5367956 | MR1 |
| -0.3675265 | -1.1455921 | -0.6310528 | MR1 |
| -0.6380936 | -0.8158955 | -0.5318305 | MR3 |
| -0.0258487 | -0.9003512 | -0.3329135 | MR1 |
| 0.7264414 | 1.0693981 | 1.0644064 | MR2 |
| -0.7606154 | -0.3005669 | -0.3748850 | MR2 |
| -0.4400121 | -0.4786269 | 0.4823075 | MR3 |
| -0.1729538 | 0.5786761 | 0.3629376 | MR2 |
| -0.4660142 | -0.0365701 | -0.5028913 | MR2 |
| -0.5160573 | -0.4888614 | -0.8377510 | MR2 |
| -0.9855923 | -0.1010805 | -0.4455918 | MR2 |
| -1.2089917 | 0.1501719 | -0.1985444 | MR2 |
| -0.2379327 | 0.3666685 | 0.5072800 | MR3 |
| -0.9078266 | 1.6250300 | 0.2234361 | MR2 |
| -0.9220901 | 0.7274965 | -0.2158877 | MR2 |
| -1.1622422 | -0.0250411 | -0.0207024 | MR3 |
| -0.9358417 | -0.6935860 | 0.0371761 | MR3 |
| -0.5404355 | 2.8733288 | -0.6636452 | MR2 |
| -1.1622268 | 0.6595541 | 0.0016911 | MR2 |
| -0.4358924 | 0.0784473 | 0.0284913 | MR2 |
| 0.1829756 | 1.8682667 | 0.6857332 | MR2 |
| -1.1591482 | -0.1504951 | -0.1140352 | MR3 |
| -0.3122178 | -0.5187054 | -0.0364263 | MR3 |
| -0.6950283 | -0.6931329 | -0.0054494 | MR3 |
| -0.7294882 | -0.3355858 | -0.5896387 | MR2 |
| -1.0263443 | 1.6731754 | -0.4002498 | MR2 |
| -1.1159914 | 1.5794419 | -0.5821503 | MR2 |
| -0.7969888 | 2.7438234 | 0.4923605 | MR2 |
| -0.4603514 | 2.9857172 | -0.7435115 | MR2 |
| -1.1496794 | 0.2091391 | 3.1950807 | MR3 |
| -1.1020217 | -0.5299483 | 0.3505327 | MR3 |
| -1.0787680 | -0.5324799 | -0.6529322 | MR2 |
| -0.3188571 | -0.5090320 | 0.0503665 | MR3 |
| -0.1753395 | 0.7040507 | 0.5297579 | MR2 |
| -0.4553281 | 0.4898921 | -0.2145291 | MR2 |
| -0.9551630 | -0.2947004 | 1.1789053 | MR3 |
| -1.0959931 | 1.3922255 | 0.1127095 | MR2 |
| -0.4422742 | 3.0019059 | -0.6339362 | MR2 |
| -2.2052603 | 8.5443991 | 3.9161805 | MR2 |
| -1.2366125 | 1.1370101 | -0.3176570 | MR2 |
| -3.3217686 | 10.5038853 | 3.9998051 | MR2 |
| -0.9847662 | 4.1832479 | 1.2487685 | MR2 |
| -0.9388075 | -1.3227448 | -0.0110649 | MR3 |
| -1.1021378 | 4.7744275 | -0.5463936 | MR2 |
| -1.0903524 | 0.2328385 | -0.4561312 | MR2 |
| -0.1180525 | 0.2188142 | -0.9213938 | MR2 |
| -0.5436459 | 1.5700404 | -0.8357402 | MR2 |
| -1.0469104 | -0.3222221 | 0.0715230 | MR3 |
| -0.1378527 | 1.4796639 | -0.4417236 | MR2 |
| -0.7265324 | 0.2440331 | -0.8507889 | MR2 |
| -0.1381405 | 2.3639121 | -0.1776378 | MR2 |
| -0.7221379 | 0.3015281 | -0.3558042 | MR2 |
| -1.0195266 | 0.5552086 | -0.4259692 | MR2 |
| -0.7945591 | 0.4599741 | -0.4984870 | MR2 |
| -0.5166759 | -0.5485588 | -0.3148669 | MR3 |
| 0.2224848 | 1.0586927 | -0.2829445 | MR2 |
| -0.6601676 | 0.8383087 | 0.6766019 | MR2 |
| -1.0070172 | 0.5229743 | -0.0404542 | MR2 |
| -1.1502283 | 0.8196942 | -0.1828275 | MR2 |
| -0.3316134 | 0.1354834 | -0.0072292 | MR2 |
| 0.3491894 | 2.3810778 | -0.2840717 | MR2 |
| -0.7943375 | -0.9834609 | 1.2523271 | MR3 |
| -0.1853149 | 0.1986097 | -0.0588846 | MR2 |
| -1.1487299 | 0.4740112 | -0.4232972 | MR2 |
| -1.2256098 | 2.1500543 | -0.1168082 | MR2 |
| -0.8365817 | 0.4647807 | -0.6504027 | MR2 |
| -0.9153940 | 2.3079543 | 0.0905829 | MR2 |
| -1.1734118 | -0.0496448 | 0.4092620 | MR3 |
| -1.2528497 | 0.5867340 | -0.7013062 | MR2 |
| -1.1146435 | -0.2568470 | 2.8068390 | MR3 |
| -0.5620033 | 0.5940168 | -1.0543719 | MR2 |
| -0.7692963 | 2.7670330 | 0.4946943 | MR2 |
| -1.1852102 | 2.0744447 | 0.2338451 | MR2 |
| -0.7193087 | 0.5267360 | 1.0571569 | MR3 |
| -0.2681011 | 1.7708039 | -0.9210606 | MR2 |
| -1.2612487 | -0.8227920 | 3.5516037 | MR3 |
| -0.6825450 | 0.4192213 | 0.0024144 | MR2 |
| -0.6540195 | 0.0779745 | 0.7048850 | MR3 |
| -1.0799796 | 0.8747948 | -0.6351929 | MR2 |
| -1.1892840 | 2.9158654 | -0.4122505 | MR2 |
| -0.9432902 | 0.9507674 | 2.1355484 | MR3 |
| -0.7783719 | -0.3451348 | -0.4765812 | MR2 |
| -0.9224370 | -0.1014693 | 0.3728341 | MR3 |
| -0.3121997 | 0.0849950 | -0.0738019 | MR2 |
| -0.2985321 | 1.5950594 | -0.4638251 | MR2 |
| 1.3806769 | 6.8602496 | -1.1114227 | MR2 |
| -0.7022650 | 0.0999142 | 1.9090310 | MR3 |
| -1.0393944 | 1.3961894 | -0.6062689 | MR2 |
| -0.1136260 | 1.2081710 | 0.2909533 | MR2 |
| -0.6802390 | -0.4543148 | -0.1593873 | MR3 |
| -0.6941744 | 0.4157125 | -0.1879515 | MR2 |
| -1.0987098 | -0.4407950 | 0.1144636 | MR3 |
| -0.0832778 | 0.9887124 | -0.6734708 | MR2 |
| -0.5780380 | -0.4878207 | -0.5810847 | MR2 |
| -0.2781742 | -0.3203201 | 1.3400103 | MR3 |
| 0.5388933 | -0.4294215 | -0.5782436 | MR1 |
| -1.2606089 | -0.2692241 | -0.6257414 | MR2 |
| -0.9369646 | -0.1244716 | 1.3305741 | MR3 |
| 0.0426521 | -0.0949825 | -0.5829014 | MR1 |
| 1.6546003 | 1.5852446 | 0.1206034 | MR1 |
| -0.0230435 | -0.4060361 | -0.8721739 | MR1 |
| -0.7659360 | -0.5693170 | 0.0978094 | MR3 |
| 1.4006359 | -1.1409185 | -0.7846840 | MR1 |
fa_pats <- loadings %>%
select(-Max, -Variable) %>%
mutate_all(as.numeric)
fa_pats <- fa_pats %>% select(sort(colnames(.))) %>% as.matrix()
print_patterns(fa_pats, pat_type = "factor", n = 1:3, title = "FA factors")## Warning: Using `size` aesthetic for lines was deprecated in ggplot2 3.4.0.
## ℹ Please use `linewidth` instead.
## This warning is displayed once every 8 hours.
## Call `lifecycle::last_lifecycle_warnings()` to see where this warning was
## generated.
## Call:corr.test(x = .)
## Correlation matrix
## MR1 MR2 MR3
## MR1 1 0 0
## MR2 0 1 0
## MR3 0 0 1
## Sample Size
## [1] 1535
## Probability values (Entries above the diagonal are adjusted for multiple tests.)
## MR1 MR2 MR3
## MR1 0.00 1 1
## MR2 0.99 0 1
## MR3 1.00 1 0
##
## Confidence intervals based upon normal theory. To get bootstrapped values, try cor.ci
## raw.lower raw.r raw.upper raw.p lower.adj upper.adj
## MR1-MR2 -0.05 0 0.05 0.99 -0.06 0.06
## MR1-MR3 -0.05 0 0.05 1.00 -0.05 0.05
## MR2-MR3 -0.05 0 0.05 1.00 -0.06 0.06
loadings %>%
ggplot(aes(x = MR1, y = MR2)) +
geom_point() + geom_label_repel(aes(label = Variable),
box.padding = 0.35,
point.padding = 0.5,
segment.color = 'grey50') +
theme(legend.position = "bottom") +
labs(title = "Variable Loadings on First and Second Factors")if ("MR3" %in% colnames(loadings)) {
loadings %>%
ggplot(aes(x = MR1, y = MR3)) +
geom_point() + geom_label_repel(aes(label = Variable),
box.padding = 0.35,
point.padding = 0.5,
segment.color = 'grey50') +
theme(legend.position = "none") +
labs(title = "Variable Loadings on First and Third Factors")
}plot_loadings <- loadings %>% select(-Max) %>% gather(key = "Factor", value = "Loading", -Variable) %>% mutate(Factor = str_replace(Factor, "MR", "Factor "))
plot_loadings %>%
ggplot(aes(x = Factor, y = Loading, fill = Factor)) + geom_col(position = "dodge") +
facet_wrap(~ Variable) +
theme(legend.position = "none", axis.text.x = element_text(angle = 45, hjust = 1)) +
coord_flip() + geom_hline(yintercept = 0, size = 0.2) +
labs(title = "Variable Loadings on All Factors")scores %>% ggplot(aes(x = Max, fill = Max)) + geom_bar() +
labs(x = "Factors", y = "Number of Individuals", title = "Number with Highest Scores per Factor") +
theme(legend.position = "none")## # A tibble: 3 × 2
## Max `n()`
## <chr> <int>
## 1 MR1 522
## 2 MR2 531
## 3 MR3 482
scores %>% gather(key = "factor", value = "score", -Max) %>% dplyr::select(-Max) %>%
ggplot(aes(x = score)) + geom_density() + facet_grid(factor~.) # should be MR1:MR4library(NMF)
num_cores <- ceiling(parallel::detectCores() / 2)
nmf_mat <- rrmc_out$L
nmf_mat[nmf_mat < 0] <- 0
res <- nmf(rrmc_out$L, rank = 3, method = "offset", nrun = 30, seed = 0, .opt = paste0('vp', num_cores))
# only compute the scores
s <- featureScore(res)
summary(s)
# compute the scores and characterize each metagene
s <- extractFeatures(res)
str(s)
W <- basis(res) # basis matrix / metagenes / contribution matrix
H <- coef(res) # mixture coeffecient matrix / metagene expression profiles / profile matrix ie. loadings ?
heatmaply(W, main = "NMF basis / contribution matrix (scores)", Rowv = F, Colv = F,
ylab = params$rowvar_name, labRow = as.character(rowlabs),
cexRow = 100, row_side_colors = data.frame("cohort" = cohorts), showticklabels = c(T, F))
heatmaply(H, main = "NMF profile matrix (loadings)", Rowv = F, Colv = F,
col_side_colors = data.frame("exposure family" = as.factor(params$colgroupings)))
loadings <- t(H) %>% as.tibble()
print_patterns(loadings, colgroups = cng, title = "NMF loadings")